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SUMMARY 


The linearized integral equation for the Foucault test of a solid 
mirror is solved by various methods: power series, Fourier series, 
collocation, iteration, and inversion integral. The case of the Casse- 
grain mirror is solved by a particular power series method, collocation, 
and inversion integral. The inversion integral method appears to be the 
best overall method for both the solid and Cassegrain mirrors. Certain 
particular types of power series and Fourier series are satisfactory for 
the Cassegrain mirror. Numerical integration of the nonlinear equation 
for selected surface imperfections shows that results start to deviate 
from those given by the linearized equation at a surface deviation of 
about three percent of the wavelength of light. Several possible pro- 
cedures for calibrating and scaling the input data for the integral 
equation are described. 
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I - INTRODUCTION 


The Foucault knife-edge test method is a widely used method of 
testing astronomical mirrors and other high-quality optical systems of 
large aperture during the final figuring. It would appear to be espe- 
cially suited for testing the mirror of a large telescope in orbit 
because of its sensitivity and minimal equipment requirements and because 
its essential requirement of a perfectly uniform test beam can be satis- 
fied in space by pointing the telescope at a star. Even if the telescope 
is free of imperfections when first built, temperature variations and 
structural relaxation or creep may subsequently introduce distortions 
so that some such method for occasionally testing the telescope optics 
is desirable. If the primary mirror is adjustable, it would be espe- 
cially desirable that the test method be capable of providing quantita- 
tive information on the distortion that must be corrected. 

This problem of the quantitative determination of the distortion 
of the mirror by means of the Foucault knife-edge test is considered 
here. The analysis is based on the work of Linfoot (Ref. l), who derived 
expressions for the distribution of the light intensity seen in the 
kinfe-edge test for both perfect and imperfect mirrors as a function of 
the mirror surface distortion. For a uniformly reflecting mirror with 
small surface imperfections, the change in the light distribution from 
that of a perfect mirror is given by a linear integral expression involv- 
ing the phase distortion of the converging wave front at the mirror. 

The solution for this phase distortion, which is directly related to the 
mirror surface distortion, involves the inversion of this integral 
expression. 

Thus, the problem of the determination of the mirror imperfections 
from the knife-edge observations reduces to the problem of solving an 
integral equation. In the work herein reported, various methods of 
solution were studied. Five function methods involving power series and 
Fourier series, a collocation solution, two iteration solutions, and an 
inversion integral solution are presented and compared for selected 
examples. Three of the more suitable methods were applied to the Casse- 
grain type of mirror, or mirror with a central hole. In addition to 
describing and discussing the advantages and disadvantages of these 
various methods, this report also presents same nonlinearized calcula- 
tions of the knife-edge test intensities for distortions of increasing 
size, in order to help the user estimate the limitations of the linearized 
solutions . 
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II - DERIVATION OF INTEGRAL EQUATIONS 


The optical arrangement in the knife-edge test is shown in the 
sketch, Fig. 2-1. The knife-edge passes through the focus of the mirror 
and the lens focuses the mirror surface onto the image plane of the lens. 
This image plane may be the location of the observers eye or may be the 
location of a group of light sensors. Linfoot (Ref. l) showed that when 
the observer looks past the knife-edge toward the mirror, the observed 
distribution of light intensity at the mirror along any line normal to 
the knife-edge depends only on the mirror figure along that line. Thus, 
with the knife-edge normal to the plane of the sketch in Fig. 2-1, only 
a line of sensors in, or parallel to, the plane of the sketch is needed 
in the image plane of the lens. To completely test the mirror, several 
lines of sensors would be needed. Katzoff (Ref. 2) discusses possible 
arrangements of the sensors for a complete test of the mirror. He also 
examines the problem of precise location of the knife-edge when the 
mirror has imperfections. 



Fig. 2-1. Optical Arrangement in Knife-Edge Test 


N. 



Fig. 2-2. Definition of x for Lines on the Mirror 
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Since the light intensity along a line across the mirror normal to 
the knife-edge depends only on the shape of the mirror along that line, 
let x he the distance along this line. Let x = -1 and x = +1 represent 
the points where the line crosses the edges of the mirror, regardless of 
whether the line is a diameter or a chord of the mirror (see Fig. 2-2). 
Linfoot (Ref. l) derived the following basic equation for D(x), the 

2jtD(x) = rtE(x) + i f — diL dt (2-1 ) 

t - X 

complex displacement at the image plane of the lens, where E(x) is the 
complex displacement of the converging wave front just as it leaves the 
mirror, and t is a running variable corresponding to the variable x. 

The observed distribution of light intensity along the line is 


I (x) = 1 k 2 Jd(x)| 2 . 


( 2 - 2 ) 


l(x) and D(x) are functions of x in the sense that the observer looking 
past the knife-edge sees l(x) as the apparent brightness at point x. 


The conditions or restrictions on Eqs. (2-1) and (2-2) are (see 
p. 138, Ref. l): 

(a) the diameter of the mirror subtends only a small angle 
at the focal point; 

(b) the errors of figure of the mirror, though they may 
amount to many wavelengths, are small compared to the 
focal distance; 

(c) the errors of slope in the mirror surface are small; 
and 

(d) the function E(x) is continuous and differentiable except 
at the edge of the mirror. 

Under these conditions, the equations are valid for mirrors of arbitrary 
edge contour, including central piercings as for a Cassegrain arrangement, 
and of variable reflecting power. 

The function E(x) can be taken in the form 


E(x) = |e(x) I exp ["- cp(x)] 

A J 


(2-3) 


where exp <p( x )J is a phase function which describes the imperfec- 


tions on the surface of the mirror; A is the wavelength of light. For 
+he perfect mirror, cp(x) = 0 and E(x) = 1, so that 
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2itD Q (x) = it + i log ( - " x ) , perfect 

\ ! + x / 


mirror . 


(2-4) 


and 

I 0 = 4it 2 |D 0 (x)| 2 = it 2 + -l<x<l, (2-5) 

where the subscript o indicates the perfect mirror. For the case of the 
distorted mirror, take 


il(x) = - cp(x)/A 

E(x) = exp[2itiri(x)3 = cos 2iti}(x) + i sin 2jtr)(x) (2-6) 


where T](x) is the phase error, in wavelengths, on the converging wave 
front where it leaves the mirror. It is twice the error, in wavelengths 
of the mirror surface at that location. The phase error ri(x) is positive 
when the imperfection on the mirror surface is raised toward the observer. 

If Eq. (2-6) for E(x) is substituted into Eqs. (2-1) and (2-2), then 


and 


2«d(x) = it [cos 2 jui(x) + i sin 2itq(x)] 


+ i f cos 2irri(t)dt f si 

J-i t * x ~J-i 


sin 2itri(t)dt 
t - x 


l(x) = it 2 + 


T 


cos 2^Ti(t)dt 
t - x 


/> 


sin 2itT)(t)dt 
t - x 


+ 2it sin 2rtt|(x) 


J cos 2n 
i t - 


2nn(t) 


dt 


r 1 

2x cos 2itn(x) I 5S_2iaMat 

JL 1 t - x 


(2-7) 


( 2 - 8 ) 


If t] (x) is sufficiently small so that 


cos 2itq(x) ~ 1.0, sin 2iaj(x) ~ 2itrj(x) 
[ rj (x) ] 2 ~ 0.0 
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then Eq. (2-8) can be linearized approximately to give 

/ i 

at (2-9) 

t - x 

This is the basic linear integral equation for tj(x). Take 

f(x) = ltot 2 t)(x), F(x) = I - I Q (2-10) 

so that the basic Eq. (2-9) is 

F(x) = f (x) log ( ) -f ■£&!. dt . (2-11) 

\ 1 + x / h t “ x 

For some of the solutions to be given later in this report, it is 
preferable to write Eq. (2-11 ) in the form 

F(x) = - f f l t ) ~ f ( x ) at (2-12) 

/ t - x 


When the mirror has a central hole of radius R, the linearization 
of Eq. (2-8) gives 


F(x) = I 


I 0 = f (x) log 


1 - x x + R 
1 + x x - R 


where 


f~ f (t)dt _ f f (t)dt 
J t-x t-x 

I 0 = Tt 2 + log 2 I 

1 1 + x x - r| 


(2-13) 


(2-14) 


for any diameter. However, as in Fig. 2-2, the line on the mirror may 
be any chord so that R in Eq. (2-13) may vary from R = 0, the case of 
Eq. (2-11) when the chord line does not cross the hole, up to R = radius 
of the hole. Thus, in Eq. (2-13 )j x = -1 and x = 1 refer to the edges 
of the mirror for any line and x = ±R, with R properly scaled to the x 
length, refer to the edges of the hole if the line crosses the hole 
(see Fig. 2-3). 
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Fig. 2-3. Knife-Edge Test Lines for Cassegrain 
Mirror System 


For later use, Eq. (2-13) can be written in the form 



( 2 - 15 ) 


where L 2 consists of the two line segments, -1.0 to -R and R to 1.0. 


The problem now is reduced to solving these integral equations for 
f(x) or q(x). Since little can be done about solving the nonlinear 
integral Eq. ( 2 - 8 ) for r|(x), the work of this report is concentrated on 
solving the linear integral Eqs. (2-ll) or (2-12) and (2-13) or (2-15) 
for f(x). Some results for the exact Eq. ( 2 - 8 ), obtained by numerical 
integration, are given in Sec. XI for assumed values of t](x). Before 
presenting the various methods of solving the linear equations, some 
discussion of the general restraint conditions, homogeneous solutions, 
numerical integration difficulties and computer limitations for these 
equations is given in the following section. 


Ill - RESTRAINT CONDITIONS ON SOLUTION OP TEE LINEAR EQUATIONS 


There are several factors that must he considered in attempting to 
solve the integral equations (2-11), (2-12), (2-13), and (2-15). 

The problems of (a) point values of F(x), (b) numerical integration, 

(c) function methods ‘with numerical integration, (d) discontinuities at 
the edges of the mirror, (e) homogeneous solutions, (f) rigid body- 
rotations of the mirror and (g) computer limitations are described below. 

(a) Point values of F(x) . Since F(x) = l(x) - I Q (x), the change 

in light intensity due to a distortion in the mirror and must be obtained 
by sensors from the knife-edge test, it mil be known only at selected 
points on selected lines on the mirror (unless it is measured with a 
scanning sensor). Thus, with F(x) given at points only, the equations 
must be solved numerically. Regardless of what method of solution is 
used, numerical integrations mil have to be made, either directly or 
indirectly. Since these equations may be singular, with the integrand 
becoming infinite when the integration variable t crosses x, there may 
be difficulties with the numerical integration. As the particular type 
of numerical integration mil depend upon the method of solution, the 
specific problems with the integration are discussed later, where the 
different methods of solution are described. Some general comments on 
numerical integration are given in item (b)„ 

(b) Numerical integration . In the numerical integration for the 
various methods of solution described later in this report, it wil l be 
assumed that readings of P(x) from the knife-edge test mil be at equally 
spaced points in a line across the mirror. For this case of equal 
intervals, the Newton-Cotes quadrature formulas (Ref. 3 5 p. 397) are 
probably the most suitable for the numerical integration. These formulas 
give weight coefficients Hj, or the Cotes numbers, to multiply the 
ordinates of the integrand for the integration. Thus 

1 J 

f (x)dx = £ Hjf(xj) (3-1) 

The two-point formula is the same as the trapezoid rule, the three-point 
formula is Simpson's rule. The Hj numbers up to a 21 point rule are 
given on pages 536-538 in Ref. 3„ The error in the numerical integration 
using these formulas is proportional to a certain derivative of the 
integrand at some point. For example, the error in the trapezoid rule 
is proportional to the second derivative; in Simpson’s rule, to the 
fourth derivative; in the nine -point rule, to the tenth derivative. 

As pointed out by Kopal, (Ref. 3), whether a higher-order mle is better 
than a, lower-order rule depends upon the behavior of the derivatives. 

This factor is important in the solution of the subject integral equa- 
tions. In some of the methods discussed later, it will be shorn that 
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the higher order derivatives involved in the numerical integration errors 
are extremely large and may he infinite at some points. This results 
in the trapezoid rule giving the best results in some methods of solution. 

(c) Function methods with numerical integration . The function 
method of solution assumes that the unknown function f (x) can be approx- 
imated by a finite sum of known functions with unknown multiplying 
constants, as 


N 

f(x) = £ a* Pi (x) (3-2) 

i=l 

The Pi (x) functions may be x 1 terms, trigometric terms or any selected 
group of functions. The constants a* are to be determined for the 
solution. The determination of these constants normally involves the 
evaluation of certain integrals containing F(x) and Pi(x). In the 
present case, these integrations must be performed numerically. This 
introduces a second approximation into the function solution in Eq. (3-2). 
Hot only is there a finite number of the a £ , they are also approximate 
and no better than the numerical integration used to evaluate them. It 
is evident that if the numerical integration is quite accurate, then 
the accuracy in f(x) in Eq. (3-2) 17111 be limited only by the number 
and type of selected functions p £ (x) . Since the larger the number of 
points used in the integration the better the results are, as many points 
as possible should be used. The number of points used depends not only 
on the behavior of f(x) but also on the selected Pi(x). If the p A (x) 
are cyclic, more points will be needed if the number of terms in Eq. 

(3-2) is large, (p. 4o8, Ref. 3). In fact, the number of points in the 
numerical integration should be at least equal to the number of constants 
in Eq. (3-2). Preferably, the number of points should be many times the 
number of constants for cyclic p^x), but from practical limitations of 
the computer, twice as many points as constants seems to be a good com- 
promise. In order to compare the various methods discussed later, the 
number of points on the mirror for all numerical integrations has been 
selected as 


40 equal intervals, 4l points for 

mirror without hole (3-3) 

40 equal intervals, 20 on each side, 

42 points for mirror with hole (3-4) 


The number of constants used in the function methods of solution has 
been selected as 


20 constants for function methods. 


(3-5) 
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Further discussion of the numerical integration in the function methods 
will be given for each of the particular function methods. 

(d) Discontinuities at edges of mirror . In Eqs. (2-11) and (2-13) 
the log terms become infinite at the edges of the mirror and the edges 
of the central hole. This may mean that f(x) must be zero at the edges 
or F(x) must be infinite at the edges. Since F(x) is a change in light 
intensity, it is not infinite at the edges. However, since the mirror 
may have a distortion at the edge, it is not desirable to make f(x) = 0 
at the edges. Thus, these factors must be considered in the solutions 
so as not to restrict the solutions more than necessary. It should be 
noted that this problem can be avoided in those methods which make use 
of Eqs. (2-12) and (2-15). 

(e) Homogeneous solutions . It is evident from Eqs. (2-12) and 
( 2 - 15 ) that 

f(x) = Ci (3-6) 

with Ci a constant, makes F(x) = 0 so that this f(x) is a homogeneous 
solution of the equations. This solution represents a uniform change in 
phase of the entire wave front and has no physical significance as the 
reference phase is arbitrary. However, the method of solution must 
include a reference point or condition for this constant. 

It can be verified that Eq. (2-15) for the hole case has another 
homogeneous solution 

*(*) - , (3-7) 

-A. 

which makes F(x) = 0. This is a possible physical distortion of the 
Cassegrain type mirror that the knife-edge test is unable to detect, as 
no change in light intensity is produced by this distortion. The evalu- 
ation of C 2 by using restraint conditions on the mirror is considered 
in Section X. 

(f) Rigid body rotation of mirror . In Eq. (2-12), it is easy to 
show that 


f(x) = C 3 X gives F(x) = - 2C 3 , (3-8) 

which indicates that a uniform change in the light intensity implies a 
rotation of the line on the mirror being observed. However , since the 
reference scale for the light intensity may be arbitrary, the rotation 
in Eq. (3-8) may be an apparent rigid body rotation of the entire mirror 
and not a real distortion of the mirror. If a group of parallel lines 
on the mirror all gave the same rotation, then it would undoubtedly rep- 
resent an apparent rigid body rotation. It is evident that the method 
of solution for f (x) in the equations should not include restraints 
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at the edges or elsewhere that would keep the method from giving the 
solution in Eq. (3-8). That is, the solution must handle the problem 
described in item (d) above. Further discussion of the scale factor 
problems on F(x) is given in Section XII. 

(g) Computer limitations . Since most of the methods studied in 
this report involve the inversion of matrices of large order (either 
4l by 4l, 42 by 42, or 20 by 20 from Eqs. (3-3) to (3-5 ) ), the question 
of the accuracy of the inversion on the digital computer must be consid- 
ered. On the IBM 360 system used, single precision arithmetic uses 
8 digits and double precision arithmetic uses 16 digits. It was found 
that the 8-digit arithmetic was insufficient for several methods, so 
that it was decided to do all the calculations in the report using 
l6-place arithmetic. Even so, one method of solution could not be com- 
pleted because 16 digits was insufficient,, 

The above listed factors influencing the solution of the equations 
w il l, be considered in each method of solution in the following sections. 
Several function methods using power series and Fourier series, a collo- 
cation method, an iteration procedure, and an inversion integral solution 
are described. In each case, the solution is set up in matrix form using 
forty (40) equal intervals across the mirror. In the solution form 

[f] = [G][F] (3-9) 

where [f ] is a 4l by 1 matrix for the mirror distortion at 4l points on 
the mirror, including edge points, and [F] is 4l by 1 for the same 
4l points (as observed in the knife-edge test on the line). The [G] 
matrix (4l by 4l) is derived as a product of various matrices for each 
method of solution. The results given by all the methods for selected 
F(x) functions as well as the [G] matrices are compared in Section IX 
for the case of no hole in the mirror. The Cassegrain type of mirror 
is considered in Section X. 


IV - POWER SERIES SOLUTIONS 


The integration in the integral equation (2-12), or 

F(x) = - f f (*> - at , (4-1) 

I t - X 

-'Ll 

is simple to perform if f (x) is assumed as a power series 

f(x) = £ A JX > (4-2) 

j=l 
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II 


whence 


F(x) = 


tJ - x^ 


t - x 


dt 


- £ a J 1 

0=1 J - i 

J f 1 

■ - £ “I “‘- 1 

0=1 L j 0 - 1 

_ V , f [1 - (-l) J - k+1 ]x k ” :L 

“ k 3 k 3 -* + 1 


f x 3-l f fl - ( -1 ) k +1 JA,, 

k k *- 3+1 


+ t J “ e x + • • • + tx J " 2 + x J " J j dt 


+ •• • + 2x s ~ 1 


(4-3) 

(4-4) 


Note that the form for f(x) in Eq, (4-2) makes the Cj. constant in 
the homogeneous solution, Eq. (3-6), zero. 

Since Eq. (4-3) has no zeros in the denominator, it should be a 
simple form to solve for the Aj constants. Three different methods of 
solution were made, as described below. 

Method A. Expand F(x) in a power series 


J 

F(x) = £ Bj-ix - 3 " 1 
0=1 

so that Eq. (4-4), for like powers of x, gives 


Bj-x = - 



(-D^ iU A k 
k - a + i 


(4-5) 


(4-6) 


The Aj constants can be obtained from the Bj constants in Eq. (4-6) by 
a back recursion, in which 
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At = - \ Bj-i 
At -i = - | St -2 
At-2 = - | [bj-s + | At] 
At-s = - | [bj-4 + | At-i] 


Aj-i+i 


1 

2 


[ 


Bj-i 


+ 


i-1 


E 


k=l 


1 - (-l) 1+1 " k 
i + 1 - k 


Aj-k+i 


i 


2,3, 


*, J 


In matrix form 


[C ] [A] = [B] 

Ctj = 0, J < i 

Ci j ■ - Q- - jW J : 1+1 i 
d - i + 1 

[A] = [erHB] 


J > 1 


To calculate the Bj constants multiply Eq. (4-5) by x i_1 
grate to get 


or 



F(x)dx 





-X 


1 - (,l)l+J-3- 
i + d - 1 


Mtl 

£ x^ 1 F(xj) HjAx = 

j=l 



Bj -l 


1 - (-l) 1 ^- 1 
i + 3 - 1 


where Hj are the Cotes numbers for the numerical integration, 
form, Eq. (4-10) is 

% [X] [H] [F] = [D][B], M equal intervals, 

M 


(4-7) 

(4-8) 
and inte- 

(4-9) 

(4-io) 

In matrix 
(4-11) 
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where [B] is a J by 1 matrix for the B c , Bj_i constants, [D] is 

J by J with the elements 


Dlj . x - , (4-lla) 

X + j - 1 

[F] is (M + l) by 1 for the value of F(x) at the (M + l) equally spaced 
points, [H] is a (M + l) by (M + l) diagonal matrix with the Hj numbers 
on the diagonal, [X] is a J by (M + 1) matrix with elements 


Xij = xj" 1 (4-llb ) 

The matrix [D] is ill-conditioned and could not be inverted for the case 
of twenty constants (j=20) using sixteen place arithmetic on the computer. 
However, the [D] matrix is similar to the Hilbert matrix. By changing 
the x variable to make the interval of integration from 0 to 1 rather 
than -1 to 1, the resulting [D] matrix is exactly the Hilbert matrix with 


Dij = 


1 

i + j - 1 


* 

From page 139 °f Ref. 4, the elements of the inverse of the Hilbert 
matrix are (matrix size J by j) 


I = (-l) i+<> (j+i-l)!(j+j -l); 

iS (i+j-l) [ (i-l) I ( j-1) ! ] 2 (J-i ) ! (J-o ) ! 


(4-12) 


For the case of J = 20, the numerator and denominator for some elements 
in Eq. (4-12) were too large for the computer to handle, so that a 
recursion multiplying factor had to be used, starting from ID i;l = J 2 . 
This gave the elements, but they varied from ID 1;l = 400 to 
XDi 5 , x 5 = 3.6(10 ) 27 for the case of J = 20. With these elements in 
[D] _1 for 


[B] = [D] -1 [X][H][F] (4-13) 

it is evident that a sixteen place computer cannot produce values of Bj 
of the order of 1.0. In fact, it gave Bj of order 10 11 as might be 
expected (order 27 minus order 1 6 = order 11). 

This result shows that a numerical function cannot be expanded 
directly into a power series with large order terms by regular computers 
with a l imi ted number of digits in the calculations. Naturally, if the 
function and its derivatives are known, then the expansion is the 
Maclaurin series with the Bj being given by the derivatives at x = 0. 


13 



Since the solution could not be obtained for a sufficiently large 
number of constants using the available computers, method A has not been 
completed in this report. 

Method B. Multiply Eq. (4-3) by x 1-1 and integrate to get 

J 

£ KijAj = E u i = J 

0=1 

where 

pi 

(4-14) 

E 1=1 x 1 " 1 F(x)dx 

(4-15) 

r _ Y Cl " - (-D i+k " 1 ] 

& (j-k+l)(k+i-l) 

(4-16) 

This method expands F(x) in a power series but does not calculate the 
Bj in Eq. (4-5) directly. Instead, the Aj are calculated directly from 
Eq. (4-l4). If desired, the Bj can then be obtained from Eq. (4-6). 

In matrix form, Eq. (4-l4) is 


[K] [A] = [E] = | [X ] [H ] [F ] 
M 

(4-17) 

where Eqs. (4-15) and (4-10) have been used. Here [A] is J by 1 for 
the J unknown Aj constants, [K] is J by J with elements in Eq. (4-l6), 
and [X], [H], [F] are as for Eq. (4-11). Now Eq. (4-2) gives 

i 1 

< 

I 1 

1 — 1 

i 1 

II 

i 1 

i i 

(4-18) 

where [f] is (M+l) by 1 for the results at the (M+l) equally spaced 
points, and [Y] is (M+l) by J with the elements 

Y ls = Xi 3 

(4-19) 

Now, if Eq. (4-17) is solved for [A] and the result put in Eq. 
there results 

(4-18), 

[f] = [G][F] 

(4-20) 

[G] = | (Y ] [K] _1 [X ] [H] 
M 

(4-21) 
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where [G] is (M+l) by (M+l) and operates on the (M+l) values of F(xj ) to 
produce the (M+l) values of f(xj) at the same points. 

For the case of J = 20 constants for the Aj in Eq. (4-2), it was 
found that the largest terms in the 20 by 20 [K] _1 were of the order of 
10 13 , which is much better than the terms of order 10 27 in Method A. 
However, with 16-place calculations in the computer, there appears to 
be some round-off error in the [K]” 1 matrix so that 20 constants probably 
represent the maximum number that can be used in a l6-digit machine for 
this Method B solution. 

Results for f(x) using Eq. (4-20) for selected F(x) functions or 
point data are given in Section IX together with a discussion of the 
results. 


Method C. Take Eq. (4-3) as 


where 


J J 

F t = F(x t ) = £ AjBj(xi) = £ B 1(J Aj 

j=l j=l 


B t i = Bj(xi) 



[1 - (-l) J - k+1 ]xj[- 1 
3 - k + 1 


i = 1,2,***, M+l 


for (M+l ) points 


(4-22) 


(4-23) 


Solve for Aj directly without expanding F(x) in a power series. This 
is a collocation solution for Eq. (4-3). In matrix form 


[F] = [B] [A] 


(4-24) 


where [A] is J by 1 for the J constants A s , [B] is (M+l) by J with ele- 
ments in Eq. (4-23), and [F] is (Mfl) by 1 for the (M+l) point values 
of F. Unless the same number of points as constants are used [B] is a 
rectangular matrix and cannot be inverted to give [A]. However, if 
Eq. (4-24) is multiplied by [B]^, the transpose of [B], then [B]^[B] is 
a square matrix with positive diagonal terms that can be inverted. Thus 


[B]T[F] = [B]T[B][A] 

win Pit o p 

[A] = [bTbI-MbFIF] (4-25) 
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If this [A] is put into Eq. (4-l8) there follows 


[f] = [G][F] 


(4-26) 


[G] = [yICbTb ]- 1 ^] 1 


(4-27) 


where [Y] is defined by Eq. (4-19) and [B] by Eq. (4-23). 

This procedure of using the transpose of the matrix to obtain a 
solution is equivalent to solving Eq. (4-24) in the sense of least 
squares. See Appendix A for the proof of this statement. 

It should be pointed out that this power series Method C is the 
only method of solution presented herein that does not involve a numerical 
integration. -However, multiplication by [B]^ is equivalent to an indirect 
numerical integration using the trapezoid rule. This operation is similar 
to the multiplication by x 1 ” 1 in Method B, followed by a numerical inte- 
gration. 

For the case of J = 20 constants for the Aj , it was found that the 
largest terms in the 20 by 20 [B^B] -1 matrix were of the order of 10 12 , 
which is an improvement over the Method B inverse matrix. Results given 
by Eq. (4-26) for power series Method C are described in Section IX. 


V - FOURIER SERIES SOLUTIONS 


By making a change of variable 


x = cos 9 , t = cos cp (5-l) 

it is possible to solve the integral equations (2-11) or (2-12) by using 
Fourier series expansions. One procedure is to change the variable in 
the power series form in Eq. (4-3) and obtain the relation between the 
constants for both f(x) and F(x) in Fourier cosine series. Katzoff has 
used this method in Ref. 2. A brief description of the procedure using 
matrices is given below as Fourier series Method B. Another procedure 
is to expand f(x) in a sine series and F(x) in a cosine series, in which 
case the integrations for Eq. (2-11) can be made directly. This method 
is described below as Method A. 
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Method A. Use the change of variable in Eq. (5-l) and take 


Since 


N 

f(x) = f(cos 9) = ^ An sin (n0) 
n=l 


(5-2) 


x = -1 gives 0 = it , x = 1 gives 0 = 0 (5-3) 

the form (5-2) restrains f(x) to be zero at the ends. If f(x) at the 
ends, then f(x) must be regarded as discontinuous at the ends, which 
requires many more terms in the series to obtain convergence. Equation 
(2-ll) now becomes 

F(0) = 2 log/ ^ in 9 - ) £ A n sin (n0) 

\1 + cos 0/ n tl 

- T A n f* sin (rap) sin 9 ^ ( 5 . 4 ) 

■'0 cos 9 “ cos 9 

By using the identity 

2 sin (ncp) sin cp = cos (n-l) cp - cos (n+l) cp (5-5) 


the integral in Eq. (5-4) can be changed to the form of the Glauert 
integral (Ref. 5, p. 92-93) 



it sin (m 9 ) 
sin 0 


(5-6) 


which occurs in thin-airfoil theory and finite-wing theory of aerodynamics. 
This gives 


f n sin (ncp) sin cp dtp _ it sin (n-l)9 - sin (n+l)9 
•*0 COS cp - COS 0 “> sin 


2 sin 0 

= - jr cos (n0) 


so that 


L- 

/ . \ N N 

F(0) = 2 log[ — sin - — ) T A n sin n0 + « Y' 

\1 + «>■ 0/^ n n il 


A n cos n 0 


(5-7) 


(5-8) 
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Expand F(0) in a cosine series by multiplying Eq. (5-8) by cos m0 and 
integrating from 0 to it . Thus 


N 


A ra + 


where 


K m „ = 


^mn A n — > m — 0,1,* **, N - 1 

o 

E m = ~2 I F(0) cos m0d0 

it 

4 / * log(_5^-«_) 

it j q \ 1 + cos 0 ) 


sin n0 cos m0d0 


(5-9) 

(5-10) 

(5-11) 


■Ao — 0 • 


The integral for K ran in Eq. (5-H) is evaluated in Appendix B to give 


K„_ = 


K__ = 


K„_ = - 


0 , for m + n odd 

m 

V 1 
na-1 * 


ft (mtn)/2 

8 Y i 

J* l 


mrt 2 £ 


m = n 


(5-12) 


rt 2 (m+n) 2_ 2k - 1 

I (m-n)/2 1 

o 1 


rt 2 (m-n) / j 2k - 1 
k = 1 


for m + n even 


Since there is no A Q term in Eqs. (5-2) and (5-9) tut F(0) may have 
a constant term, it is necessary to use care in writing Eq. (5-9) in 
matrix form so as to properly include the m = 0 case. If the m = 0 
equation in (5-9) is put in the last row rather than the first row, then 
the matrix equation takes the simple form 


[I + K] [A] = [E] 


(5-13) 


where [I] is the N by N identity matrix with element Ijjpj = 0, [K] is 
N by N with the elements in Eq. (5-12) but with row N having m = 0 or 
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1 


(5-14) 


Ktf„ = 0 for n odd 
n/2 

16 


rot 2 k = i 2k - 1 


for n even 


and [E] is N by 1 for the values of the integral in Eq. (5-10) for 
1 = 1, 2,*", S - 1, 0, respectively. [A] represents the N unknown 
constants . 

The matrix form of E ra for the numerical integration for point values 
of F(0) is 


[E] = -L [C][H][F] 

it- 


(5-15) 


where [C] is N by M (N constants and M points for the numerical integra- 
tion) with elements 


Cij = - cos i0j , cos N0j = -1 , (5-16) 


[H] is the M by M weighting matrix for the numerical integration, and 
[F] is M by 1 for the M point values of F. Here 0 X = it for x = -1 so 
that the numerical integration starts at the upper limit in Eq. (5-10). 

Put Eq. (5-15) into Eq. (5-13) to get the constants A n in Eq. (5-2) 

[A] = 4 [I+K] _1 [C][H][F] (5-17) 

it 

The matrix form of Eq. (5-2) is 


[f] = (S ] [A] 

where [S ] is M by N with elements 


S u = sin J0i 


Finally, 


If] = [G][F] 


(5-18) 

(5-19) 
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( 5 - 20 ) 


[G] = -% [S][I+K]“ 1 [C][H] 

it 


where [G] is M by M and operates on the M values of F(0j ) to produce 
the M values of f(0j)» 

Equal AX and Equal A0 Intervals . Since a change of variable 
(Eq. 5-1) is used in the Fourier series analysis, there arises the ques- 
tion of what intervals to use in Eqs. (5-19) and (5-20). If equal Ax 
intervals are used, then the A 0 intervals will be unequal. On the other 
hand, if equal A® intervals are used, then the Ax intervals will be 
unequal. Since the other methods of solution presented herein use equal 
Ax intervals, it is desirable to calculate f(xj) at equal Ax intervals 
in order to make a direct comparison with the other methods in Section IX. 
This was done by using 0£ in Eq. (5-18) at the equal AX intervals. On 
the other hand, the numerical integration involving Eq. (5-15) can be 
done either way by using Qj in Eq. (5-16) for equal Ax or equal A 9 with 
the corresponding values of F (0j), either calculated in selected cases 
or interpolated in actual cases. The case of equal AX is given in 
Section IX for this Method A using the trapezoid rule for the matrix [H]. 


Special Numerical Integration Procedure . As pointed out in paragraph 
(c) of Section III difficulties may occur in the numerical integration 
for function methods, particularly for cyclic functions. For the case 
of 20 constants and 4l points used in this report, it was found that the 
trapezoid rule for the matrix [H] in Eq. (5-20) gave poor results for 
both equal Ax and equal A0.(see Sec. IX for AX results). In order to 
improve the numerical integration in Eq. (5-10), the cosine function was 
integrated across each interval (both for equal Ax and equal A©) with 
F(0) being held constant over the interval. Thus 


(0i+0 i+ i )/2 


,/ w 

( 01 + 01 - 0/2 


_ 

F(e) cos mede = F^O p 2 3 . 111 m9 J 

L m J fa. 


(0i+0i+l )/2 


(0i+0 i -x )/2 


= — ("sin £ (®i+0i+i ) - sin ~ ( 0i+0i-;i)l (5-21) 

m L 2 2 J 


2 2 
The matrix [C][H] in Eqs. (5-15) and (5-20) is replaced by a matrix 


[D] = [C][H] 


( 5 - 22 ) 


where, from Eq. (5-21), the elements of [D] are (0 X = it, 0^ = 0) 


20 



Dm = \ (rt-0 2 ) , Dnm = | 0M-1 

dnj = - — (0j+i-0j-i) ? d = 2,3,**-, m - i. 

Dim = \ sin^l 0M-i) , i = 1,2,3,- N - 1 

Du = (“l) i DiM , i = 1,2,3, ••*, N - 1 
Dij = - ^ j^sin — (9j+0j+i) - sin ^ (0j+0j-i)J> 
i/H, j / 1, M 

The [G] matrix in Eq. (5-20) becomes 


f (5-23) 


[G] = -% CS][I+K] -1 [D] (5-24) 

rt 

where the elements of [S], [K], and [D] are given in Eqs. (5-18), (5-12), 
(5-14) and (5-23), respectively. Also, 1^ = 0 in [I], 

The matrix [D] in Eq. (5-24) can be calculated for equal Ax or for 
equal Ad- Results for the ax case are shown in Section IX for Method A. 

Method B . This method uses the change of variable in Eq. (5-1) in 
the power series method after integration, Eq. (4-3), and relates the 
Fourier cosine series constants for f(x) and F(x). The details of the 
procedure is given by Katzoff in Ref. 2. A brief outline of the method 
in matrix form is given below so as to set up a form similar to the 
other methods in order that results of all the methods can be compared 
directly. Put Eq. (5-1) into Eq. (4-5) to get 

J J-l 

F(e) = Z B s~ i( cos 0 ) 3_:L s B o + £ ( cos e) J Bj (5-25) 

0=1 0=1 


Take the Fourier cosine expansion of F(0 ) as 

p J-l 

F(9) = -°+ £ Pj cos j0 
2 0=1 


( 5 - 26 ) 



where the Pj may be regarded as known. 


Pj = - [ F(0) cos j0d0 (5-27) 

It 


Now the Bj constants in Eq. (5-25) can be obtained in terms of the 
known Pj constants in Eq. (5-26) by using the summation form of Eq. (7) 
in Ref. 2: 


K 

cos n0 = 

i=0 


(-l) i n(n-i-l) ! (2) n ~ gi ~ 1 (cos 0) n " 2i 
(n-2i)!i! 


K = n -~.. 1 , for n odd , 


(5-28) 


K = - 


for n even 


Substitute Eq. (5-28) into Eq. (5-26) and equate the result for F(0) to 
the form for F(0) in Eq. (5-25). The coefficients of corresponding 
(cos 0) J terms on both sides of the resulting equation must be equal so 
that 


J-l 


B c = ~ + | L [1 - (-1) J+1 ](-1)V 2 P d 


2 A 


B m = E 


J-l (-l)(j-»)/ 2 [l - (-l)J+"+*]2 m - g (j) ( j + g ~ 2 ) ! 


Pj 


i(5-29) 


D=m 




m = 1, 2, • • • , J - 1 


In matrix form, Eq. (5-29) is 


[B] = [R][P] 


(5-30) 
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•where the elements of [R] are 


Rij = 0 , j < i , 


Ri 


s = J- Tv ’ 


\ (5-31) 


3 > i , 


Rll = i 
2 


The Aj constants in Eq. (4-2) are related to the Bj constants by 
Eq. (4-8) as 


[A] = [Q]" 1 ^] 

where the Qu elements are 

Qij = 0 , j < i 


(5-32 ) 


q„ . - a. - , j > , 

j - i + 1 


(5-33) 


Hence, the Aj constants can be expressed in terms of the Pj constants 
by putting Eq. (5-30) into (5-32) to get 


[A] = [QD-HRHP] 


(5-34) 


The numerical integration for the Pj constants in Eq. (5-27) is similar 
to that for the Ej constants in Eqs. (5-10) and (5-15) for Method A. 

Thus 

[P] = | [C][H][F] = § [D][F] (5-35) 

J*. It 

for the special form in Eqs. (5-22) and (5-23)* However, the last row 
in [D] should be transferred to the first row in the [D] for Eq. (5-35) 
in order that the [P] matrix be compatible with the [R] matrix in 
Eq. (5-34). 
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Since it is desirable to calculate f (xj ) at equal Ax intervals in 
order to compare the results to the other methods, Eq. (4-2) or (4-l8) 
can be used for f (x) to give 


[f] = [G][F] (5-36) 

[G] = £ [ Y ][q]-1[R][d] ( 5-37) 

rt 

where Eqs. (4-l8), (5-34), and (5-35) have been combined. Here, the [Y] 
matrix should be calculated from Eq. (4-19) for equal Ax intervals, but 
the [D] matrix should be calculated for equal A0 intervals (Eq. (5-23) 
with row N used as row 1. This requires that [F] in Eq. (5-36) be known 
at equal A 0 intervals. If the Fj values are given at equal Ax intervals, 
then an interpolation matrix can be used to calculate Fj at the equal 
A 0 intervals. For this case [G] in Eq. (5-37) becomes 


[G] = 2 [y][q]“ 1 [r][d][im ] (5-38) 

Tt 

with [IM] an interpolation matrix. Results are given in Section IX for 
a straight line interpolation matrix and for exact Fj at equal A 0, using 
a selected F function. 

To calculate fj at equal A0 intervals, as Katzoff does in Ref. 2, 
a further modification can be made to the [G] matrix in Eq. (5-37). 

If f (0) is expanded in a Fourier cosine series without the constant term, 
as 


J 

f(9) = £ Hj cos j0 , (5-39) 

j=l 


then a relation between the Aj constants in Eq. (4-2) and the Hj constants 
can be obtained in the same manner as above for the relation between the 
Bj and Pj constants in Eq. (5-30). Thus 

[A] = [S][H] (5-40) 

where 

S 13 


SiJ 


= 0 


for j < i 


(-l)(j-0/2 [l - (-l) 3+i+1 ]2 i - 2 (j)( 3 - + * ~ 




(5-4l) 
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Pat Eqs. (5-40) and (5-30) into Eq. (5-32) to get 

[S][H] = [q] _1 [r][p3 

or (5-42) 

[hi = [sr^Qr^RHP] 

It can be shewn that 

[S ] -1 [Q]“ 1 [R] = | [Q]" 1 (5-43) 

so 

[H] = ^ [Qj-'-tP] (5-44) 

2 

which corresponds to the relation given by Katzoff in Ref. 2. This 
Eq. (5-44) gives the constants in the Fourier cosine series for f (9) in 
terms of the constants in the Fourier cosine series for F(0). Thus, 
for fj at equal A 0 intervals 


[f t ] = [cos 1 [Hj ] (5-45) 

whence 

[f] = [G][F] 

[G] = ^ [T ] [Q]“ 1 [D] (5-46) 

where the elements of [T] are 

T il3 = cos j-0 1 (5-47) 

and [Q]" 1 and [D] are the same as in Eq. (5-37). 

In Ref. 2, Katzoff uses the same number of points as constants, 
does not include the end points, and uses the trapezoid rule so that his 
[D] matrix has the elements 


Did = cos (i-l)0j , i,j = 1,2,3, N (5-48) 


with [D] being a square matrix. 
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VI - COLLOCATION SOLUTION 


In Sections IV and V, function methods of solution of Eqs. (2-11) 
or (2-12) using power series and Fourier series have been described. 
Another method of solution, which is simple to apply, is a direct collo- 
cation on Eqs. (2-11 ) or (2-12). Take Eq. (2-12) in the form 

F( Xi ) = -[ iM ~ f ( x *) dt (6-1) 

J-! t " x i 


which gives M equations for M X* points. Now, a numerical integration 
can be made on the unknown function f(t). Since f(xi) occurs under the 
integral, it is necessary to use the same points in the numerical inte- 
gration as in the Eqs. (6-1) in order to get a determinate system (see 
p. 455, Ref. 3). 

The question arises as to which of Eqs. (2-11) or (2-12) to use in 
the collocation. It is evident that the edge points cannot be used in 
Eq. (2-11), as the log term becomes infinite at the ends. Also, the 
integrand in the integral in Eq. (2-ll) becomes infinite at t = x. 
However, this problem can be handled by using the Cauchy principal value 
for the integration. It was also found that the collocation solution 
using Eq. (2-11 ) gave a rigid body rotation to the f(x) solution (item 
(f). Section III). This rotation apparently resulted from the first term 
in Eq. (2-11) being exact, while the integral was approximate. Now, if 
f(x) is continuous, these problems do not arise in Eq. (2-12). This can 
be shown by expanding f(t) into a Taylor's series about t = x: 

f(t) = f(x) + f ' (x) (t-x) + (t-x) 2 +••• (6-2) 

2 ! 

This gives the integrand in Eq. (2-12) as 

g(t,x) = f C x ) = f'(x) + f (t-x) +••• (6-3) 

t-x 2! 


Thus, for t = x, 


g(x,x) = f ' (x) , t = x 


(6-4) 


and the integrand is finite and continuous at t = x, if f T (x) is 
continuous. 
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Further, from Eq. (6-3) 


g’(x,x) = | f"(x), g"(x,x) = | f"’ (x), t = x . (6-5) 


This indicates that for any derivative of f(x) that is discontinuous, 
the next lower derivative of g(x,x) will he discontinuous, and the cor- 
responding derivative of g(x,x) will be infinite. As pointed out in 
paragraph (b) of Section III, the behavior of these derivatives will 
affect the numerical integration in Eg.. (6-1). Further discussion of 
these derivatives in the collocation solution is given for selected 
functions in Section IX. 

On the basis of the above discussion, Eq. (2-12) or Eq. (6-1) will 
be used in the collocation solution. With numerical integration, 

Eq. (6-1) becomes 

M 

X) (ax)Hj f ( x j) ~ f ( x i) = _ F(x t ) (6-6) 

j=l X J - x i 


where the Hj are the Cotes numbers for equal intervals in the integration 
(paragraph (b). Sec. III). For the M points and M-l intervals on 
x = -1 to x = 1, 


xj = (j-l)Ax -1, Xi = (i-l)Ax -1 


Ax = 


M-l 


and Eq. (6-6) beccmes 
M 


Hj f j .~ = - Fi ; i 


j=l 3 - 1 
As j ->i, Eq. (6-4) gives 


= 1,2, * * * , M; j^i 


(6-7) 


(6-8) 


f .J — I_. l l » f » (Xi )ax 

3 - i 


= f i-l , i ^ 1,M 


(6-9) 
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At the end points, i = 1 and M, the three point Lagrangian differentiation 
formula (p. 516, Ref. 3) gives 


f ' _ ~ 3fi + jtfg ~ f 3 
1 " 2 (Ax) 

f« = f M-g_- ffifr-i +. 3fM (6-10) 

2(Ax) 


These Eqs. (6-8)-(6-10) can be combined into a matrix equation 


[A] [f ] = - [F] 


(6-H) 


where [A] is M by M for M points with the elements 
A 12 = 2 H x + H 2 , A 13 = | (Hs-Hx) , 

= -2H m - H m _x , A M?M _ 2 = i (H m -H m _ 2 ) , 

= - Hi _ x - | Hi 
Aiji+i = Hi+i + ^ Hj 

Au = -JLl — , i ^ j , except above values, 

3 -~i 
M 

An = - £ A u , see Eq. (6-8) . 

j=l,j/i 


% 


> 


( 6 - 12 ) 


Since the diagonal elements An are the negative sum of all other elements 
in each row, it is evident that when f(x) is constant, F(x) will be zero 
as required by the homogeneous solution in paragraph (c). Section III. 

This result is not given exactly by Eq. (2-11), since the log term on the 
diagonal is only approximated by the negative sum of the other elements 
in each row. 

Since all rows of the [A] matrix in Eq. (6-12) add to zero, the 
matrix [A] is a singular and cannot be inverted. It is necessary to add 
a boundary condition on f (x) to obtain the homogeneous solution in 
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El- (3-6). To correspond to the power series solution in Section IV, 
which made f (0) = 0, the restraint condition 


f(0) = f(M+l)/2 = 0 


(6-13) 


was selected. Now the system of equations in (6-11) and (6-13) bas M + 1 
equations in M unknowns. To solve such a system, multiply by the matrix 
transpose, as in power series Method C above. In this case, the resulting 
system of M equations for the M unknowns can be expressed as 


[B][f] = - [A] T [F] (6-l4) 

[B] = [A] T [A] + [C] (6-15) 

where the elements of [C] are 

Cij =0 , C(M+l)/2i( m+i)/2 = 1*6 (6-16) 

The [B] matrix is nonsingular so that 

[f] = - [B] -1 [A]^ [f] = [G][F] (6-17) 

[G] = - [B] _1 [A]^ (6-18) 


Results given by Eq. (6-17) for selected functions of F(x) are 
shown in Section IX and compared to the results of the other methods. 
The [G] ma trix in Eq. (6-l8) is also compared to the other [G] matrices 
in Section IX. 


VII - ITERATION SOLUTIONS 


Iteration is a classical procedure for solving integral equations . 
However, when numerical integration is involved in the iteration and 
the equation has singular points, difficulties can arise. As pointed 
out in Section VI on the collocation solution Eq. (2-11 ) has infinities 
in the log term at the end points and has a singularity at t = x. Yet 
Eq. (2-11) is in the form needed for iteration so that it should be used 
rather than Eq. (2-12). 
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Equation (2-ll) can be written in two forms for iteration 


f.W 


j , — r f(I) + f 1 

1 - x \ t - X 

1 + X / 


(7-1) 


t - X 


= - F(x) + f ra -i (x) log(^-f) (7-2) 


= - G«-x(x) 


The form in Eq. (7-l) requires a direct integration to get the next 
f m (x) while the form in Eq. (7-2) requires the solution of a well-known 
integral equation for which the inversion integral is known. The form 
in Eq. (7-l) was found to be divergent so that the following discussion 
is restricted to Eq. (7-2). 

The solution of Eq. (7-2) for f m (x) is given as an example in 
Section VIII. From Eq. (8-40) 


n J_ x X(t)(t-x) 


(7-3) 


where the possible forms of X(x) and $k(x) ^ rom Eqs. (8-4l) — (8-44) are 


Xx (x) = — L= , $ h (x) = C ° 

VI - x 2 1 Vl - x 2 


f (x) may be unbounded at x = ± 1; 


X 2 (x) 


/ I - x 

1 + X 


, $ho( x ) = 0 


f m (x) may be unbounded at x = -1; 


(7-4) 


(7-5) 


X 3 (x) = W ^ + * , $h 3 ^ x ) = 0 » 


f m (x) may be unbounded at x = 1; 


(7-6) 


30 


X 4 (x) = n/1 - X s , $h 4 ( x ) = 0 > / 

1 (7-7) 

f m (x) is bounded, but Eq. (8-45) I 
restricts G m _x (t ) . 


Since f (x) is actually bounded in the problem under consideration, any 
of the four forms in Eqs. (7-*0-(7-7) may te used for X(x) and X(t) in 
Eq. (7-3). For the numerical integration in Eq. (7-3) > "the form Xx(x) 
will make the integrand zero at the ends while the other forms will make 
it infinite at one or both ends. Necessarily, the end points must be 
omitted in the calculations, but the total area in the integral must be 
included. A reasonably accurate representation of the area near the ends 
can be obtained for the X x form, but only a rough approximation can be 
obtained for the other forms. The matrix form for all the cases is as 
follows . 


Start with f 0 = 0 as the first approximation in the iteration 
procedure so that 



F(t)dt 

X(t)(t-x) 


(7-8) 


where $^(x) = 0 for all cases. 


In matrix form 


[fi] = \ [XplCAilCXp]- 1 ^] (7-9) 

JT 


where the same points are used for x and t and 


X D (i,j) = 0 , X D (i,i) = X(xi) (7-10) 


A I (i,j ) = Ai j in Eq. (6-12) , i ± j 

Aj(i,i) =0 , A 12 = H 2 , Ax 3 = i H 3 

A-MsM-l = “ h m-i j A m , m _ 2 = - H m _2 


(7-11) 


The numerical integration matrix [Aj ] is taken in the same form as that 
for the collocation method in Section VI. 


31 



For the iteration from Eq. (7-2) and Eq. (7-3) 


fa(x) 


X(x) 

f F(t)dt 

r fx(t) iog( 1 + t ) M 

rt s 

Jli x(t) (t-x) 

J x X(t) (t-x) 



Mt) i°g(^ 4)dt 
X(t)(t-x) 


or 

[f 2 ] = A. [Xd]^]^]- 1 ^] - ^ [XdJCA!]^]^]' 1 ^] 
nr nr 

= [i [Xe ]^]^]- 1 - A. [X D 3 1 A X ] [L ] E A X ] [x D ] - 

= \ [X D ][I + J ] [Aj ] [X d ] -1 [F ] 

■Si 

where [i] is the identity matrix and 


[J] = - \ [Aj ] [L ] 

it 



The form for f 3 (x) is 


^3(x) = fi(x) 


X(x) 

* 2 


r 

•'-i 


f 2 (*) K«(rfT)it 


X(t) (t-x) 


or 


[fa] = \ [X D ][I + J + J 2 ][A i ][X I) ]” 1 [F ] . 

jt 

After m iterations. 


[f»] = [S,][F] 

[G„] = % [X D ] [I + J + J 2 +•••+ J m_1 ] [Aj ] CX D ] -1 

It 


( 7 - 12 ) 

[F] 

(7-13) 

(7-14) 


(7-15) 

(7-16) 
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It should be noted that as m -» » in Eq. (7-16), 


[G»] -»[G] = [X D ][I - m -»« . (7-17) 

IT 

However, this l imi t value of [G] is just the [G] for a collocation 
solution of Eq. (7-3) with = 0, f m = f ra -i. That is, Eq. (7-3) becomes 

f 1 F(t)at f 1 t.w rrl)at 

X(x) J_ ± X(x) (t-x) J_ x X(t)(t-x) 


or 


[Xur 1 ^.] + \ [A I ][L][X D ]" 1 [f m ] = \ [Aj][Xj)]“ 1 [F ] , 

it rt 

or, hy use of Eq. (7-l4), 

[i - j][x D ] _;L [f ra ] = ± [AjJCXd]- 1 ^] , 

whence 

[f ra ] = [G][F] (7-18) 

where [G] is given hy Eq. (7-17). This indicates that the iteration 
solution can only approach a collocation solution of Eq. (7-3) with 
G m _ x (t) given by Eq. ( 7 - 2 ) and f m -i = f m at convergence. Because of the 
log term and the X(t) term in Eq. (7-3) 5 it appears that a collocation 
solution of Eq. (7-3) would be less accurate than the collocation solution 
in Section VI for the original Eq. (6-1). Thus, there is little justi- 
fication for using a collocation solution on a modified equation or for 
a iteration solution in this particular problem, as compared to a collo- 
cation solution of the original equation. However, to get some idea of 
the behavior of the iterations and of the different forms of the X(x) 
function, the Eqs. (7-15) and (7-16) are used for five iterations on a 
selected function in Section IX for the X! and X 2 forms, Eqs. (7-4) and 
(7-5 ) ) and the results compared to the other methods. For the X 2 form 
in which the integrand in Eq. (7-3) becomes infinite at x = 1, the [Aj] 
matrix in Eq. (7-11) was modified to include the approximate integrated 
area over the interval adjoining x = 1. Examination of the powers of 
the [J] matrix in Eq. (7-16) shows for both Xj. and X 2 that up to the 
fourth power, most of the elements axe decreasing but some are oscillating 
and some are actually increasing. Because of the approximations in the 
numerical integration it does not appear that al l the elements in J™" 1 
will become zero for a practical number of m iterations. This error 
accumulation can be avoided by using the limit form in Eq. (7-17 ) > but 
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then this is not an iteration solution hut a collocation solution. 
Further discussion of the iteration results is given in Section IX 


VIII - INVERSION INTEGRAL SOLUTION 


It is possible to obtain a solution for both the solid mirror, 

Eq. (2-11), and the mirror with a central hole, Eq. (2-13), in an integral 
form by using methods described by Muskhelishvili in Ref. 6. Although 
the formula from which the solutions can be obtained is given on page 328 
of Ref. 6, the definitions, conditions, and derivations leading up to 
this particular formula are spread throughout the book. The derivations 
are made for the complex plane and the equations are in terms of the 
complex variables z and t. A brief outline of the derivation of the 
solution is given below. 

Consider the limiting value of the Cauchy integral 

®(s) = -k. / (8-3.) 

2jri l t - z 

as z approaches t Q on the arc or contour L. If cp(t) satisfies a Holder 
condition on L of the form 

|cp(t 2 ) - cp(t x ) | < A|t 2 - tujk , H condition, (8-2) 

where A and k are positive constants, then Muskhelishvili (Ref. 6) shows 
that $ (z) is continuous on L from the left and from the right, with the 
exception of those ends at which <p(t) ^ 0. Further, <f(z) tends uniformly 
to the limits 


„ +( t o) = | <p(t 0 

(8-3) 

.-(to) --!»(*„)♦£ 

(8-4) 


from the left of L and from the right of L, respectively. Here, L may 
be a union of smooth, non-interesecting arcs or contours, L x , L 2 -**Lq, 
with definite positive directions with the + region on the left. These 
equations can be written in the form 
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$> + (to) - «"(t 0 ) = <p(t Q ) 


(8-5) 


» + (to) + t-(to) - i- J = 2* (to) (8-6) 

** L t " *0 


The f “unction $ (z) , holomorphic in the complex plane, except possibly at 
infinity, and continuous on L from the left and from the right with the 
possible exception of the ends b k , near which the inequality 

b(z)l < J cons ~ fp ,0<P<1 (8-7) 

|z - bkP 

holds, is called sectionally holomorphic with the line of discontinuity 
L or boundary L. 

From Eq. (8-5) , it is evident that the Cauchy integral formula 
solves an important boundary value problem. That is, it is required to 
find a sectionally holomorphic function o(z) , zero at infinity, and 
satisfying the given boundary condition 


$ + (t) - $“(t) = cp(t) on L . (8-8) 

If cp(t) satisfies the H condition, Eq. (8-2), the solution for $(z) is 
given by Eq. (8-l). 

In his book on "Singular Integral Equations, "(Ref. 6) Muskhelishvili 
develops this procedure of solving boundary value problems and applies 
it to various boundary value problems and integral equations . To arrive 
at the integral equation form that is suitable for the subject problem 
(Eqs. 2.11 and 2.13), it is necessary to first consider the Hilbert 
boundary value problem for both the homogeneous and nonhomogeneous cases. 


Homogeneous Hilbert Problem 


The homogeneous Hilbert problem can be stated as : To find a 

sectionally holomorphic function $ k (z) of finite degree at infinity 
satisfying the boundary condition 


4>h(t) = G( t )®h( t ) on L ’ (8-9) 

where G(t) is a given non -vanishing function on L, satisfying the H 
condition. 
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Take logarithms on both sides of Eg. (8-9) to get 


log $£(t) - log $ k (t) = lo & G ( t ) 


( 8 - 10 ) 


which has the form of Eg. (8-8) above. Thus, the formal solution for 
log is given by Eg. (8-l) as 



log ®h(z) 

if log 0(t)dt 
2jti £ t - z 

(8-11) 

or 

$ h (z) = e p ( z ) , 

P ( z ) = 1 / lo S dt 

2ni t t - z 

(8-12) 


There are two difficulties with this solution for $ k ( z ) in Eg. (8-12). 
The function log G(t) is multivalued and the solution for may not 
obey Eg. (8-7) at the ends of the arcs Lj (j = l,***g) of L. These 
difficulties are resolved by Muskhelishvili (pp. 230-234, Ref. 6 ) in 
the following manner. 

Near any end b^(k = 1,**», 2g), P(z) in Eg. (8-12) can be written 
in the form 


P( z ) = + l0g G ^ k ' > log (z-b k ) + P Q (z) 
2iri 

where P Q (z) remains bounded near b k . Thus 

e p ( z ) = (z-b k ) ak+iCk q(z) 

near b k , where a k and c k are real constants given by 

a k + i= k = t 108 G(b ^ 

2jti 


(8-13) 


(8-14) 


(8-15) 


and Q(z) is a non-vanishing function bounded near b k . Now select 
integers A k , satisfying the conditions 


and put 


-1 < a k + A k < 1 

T(z) = (z-b 1 )^ 1 (z-b 2 )^ 2 - • *(z-b 2 g )^ 2(1 


(8-16) 

(8-17) 
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Then the function 


X(z) = T(z)e p ( z ) 


(8-18) 


satisfies Eqs. (8-7) and (8-9) and is a solution of the homogeneous 
Hilbert problem. 

Since Eq. (8-l6) may permit two values of A^ for some k values , 
the solution X(z) is not necessarily unique. Its form will depend upon 
what requirements are placed upon the solution at the end points. If 
the solution is to be bounded at. certain end points * then one form 
results. If it is permissible for the solution to become infinite at 
some end points (under condition (8-7) ) 5 then a different form of the 
solution results. Muskhelishvili (p. 231 ? Ref. 6) divides the solutions 
into classes on basis of the end point requirements and shows that the 
most general solution in a certain class has the form 


$ h (z) = X(z)R(z) 


(8-19) 


where R(z) is a arbitrary polynomial of order depending upon the desired 
behavior of the function at infinity. Thus the solution of the homogen- 
eous Hilbert problem has the form 


$h( z ) = T(z)R(z)e P ( z ) 


( 8 - 20 ) 


where P(z) and T(z) are defined by Eqs. (8-12) and (8-17)> respectively, 
and R(z) is a polynomial with arbitrary constants. 


Nonhomogeneous Hilbert Problem 


The nonhomogeneous Hilbert Problem may be stated as : To find a 

sectionally holomorphic function $>(z), having finite degree at infinity, 
for the boundary condition 


3> + (t) = G(t)<jT(t) + g(t) on L 


( 8 - 21 ) 


where G(t) and g(t) are given functions with at most finite discontinui- 
ties on L and G(t) ^ 0 on L. As for the homogeneous problem, L may be 
a union of smooth nonintersecting arcs. 

The solution of the nonhomogeneous problem can be obtained from the 
homogeneous solution. From Eq. (8-9) 
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„ m . 4W = Slot) 

<&£(t) X-(t) 


( 8 - 22 ) 


where X(t) is $ h (t) in Eq. (8-20) with the unknown polynomial R(t) 
omitted, or 

X(t Q ) = T(t 0 )e P ( t o) (8-23) 

With cp = log G and $ = log X, it follows from Eqs. (8-5) and (8-6) that 


X + (t D ) = G(t 0 ) X(t D ) , 

X'(t 0 ) = X (t Q ) / n/ G(t 0 ) . 


Put G(t) from Eq. (8-22) into Eq. (8-21 ) to get 

$ + (t) _ $"(t) _ g(t) on L 
x+(t) x-(t) x*-(t) 


(8-24) 


(8-25) 


This Eq. (8-25) is of the form of Eq. (8-8) so that from Eq. (8-l) 


$ 


(z) = _ 1 _ [ 

/ \ ^ • J 


g(t) 


X(z) 2rti £ X + (t ) (t-z) 


dt 


Thus 




(8-26) 

(8-27) 


where any homogeneous $h( z ) has been included. By using Eq. (8-24), 
the solution (8-27) can be written as 


$(z) = f dt + $h(z) 

2rti l X(t)v/G(t)(t-z) 


(8-28) 


which is the general solution of the nonhomogeneous Hilbert problem with 
X(t) given by Eq. (8-23). 

This solution (8-28) can be used to solve various singular integral 
equations, including the integral equations (2-11) and (2-13) under 
investigation in this report. 
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Integral Equations Solutions 


The integral equation discussed in Section 107 of Ref. 6 can be 
adapted to the problems involved in this report. The equation has the 
form 


A(t 0 )f(t 0 ) + d t = F(t c ) , on L, (8-29) 

rtl . l % ~ 


with A(t), B(t), F(t) satisfying the H condition, Eq. (8-2). If f(t Q ) 
corresponds to <p(t c ), then Eqs. (8-5) and (8-6) can be used to express 
the unknown f(t D ) in terms of an unknown function $(z) and the Eq. ( 8 - 29 ) 
can be written in the form 


or 


A(to) [<$ + (to) - $”(to)] + B(to)[$ + (to) + $“(to)] 


= F(t 0 ) on L 

( 8 - 30 ) 

$ + (to) “ G(t 0 )<I> (t G ) + ^ 

A(t c ) + B(t 0 ) 

( 8 - 31 ) 

G(t 0 ) = A(t °) " 

A(t 0 ) + B(t 0 ) 

( 8 - 32 ) 


This Eq. (8-31) corresponds to the nonhomogeneous Hilbert problem 
in Eq. (8-21 ) so that if $(z) is determined for this problem then the 
solution of the integral equation ( 8 - 29 ) is simply 


f(t D ) = ®+(t Q ) - ®-(t Q ) 


(8-33) 


From Eqs. (8-28), (8-3l) s and (8-32) 

®(z) = { -= lCt)dt + ® h ( z ) 

2rti L ^G(t ) X(t)[A(t) + B(t)](t-z) 


= X(z) f F(t)dt + 

2iti { X(t )' n/a 2 -B 2 (t-z) 


«h(z) 


(8-34) 
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Use Eqs . (8-3) an d (8-4 ) for $(z)/x(z) in Eq. (8-34) and take t 0 = x, 
cp(t) = F(t)/X(t)*VA 2 -B 2 , to get 

4 +( x ) = 1 X + (x)E(x) + X+(x) f F(t )dt 

2 X(x)n/a 2 -B 2 2rti L X(t )V a 2 -B 2 (t-x) 

+ ®J(x) (8-35) 

$ -( x ) = . 1 X~(x)F(x) , X~(x) f F(t)dt 

2 X(x)'Va 2 -B 2 2rti L X(t)\/A 2 -B 2 (t-x) 


$h( x ) 


(8-36) 


Since, from Eq. (8-24), 


x- + x-JJg + i)x= -M= 
\ sTg) *7a^F 


X+ - X" = 


-2BX 


/a 2 - 


B" 


- K = (X + -X“)R(x) = - 


2B 


$h 


it follows from Eqs. (8-33)> (8-35) and (8-36) that 

f(x) = A(x)F(x) _ B(x)x(x) f F(t)dt 

A 2 (x) - B s (x) 3tWA 2 -B 2 L X(t)'s/A 2 -B 2_ (t-x) 


(8-37) 


- Mil- $ h (x) (8-38) 

■Va 2 -b 2 


This is the general solution of the integral equation (8-29) with $^(x) 
given by Eq. (8-20) and X(x) given by Eq. (8-23). The solution can be 
used for the particular integral equation used in the iteration method 
in Section VII as well as for the integral equations (2-ll) and (2-13). 
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Integral Equation Solution for Iteration Method 


In Eq. (8-29) take A(t 0 ) = 0, B(t 0 ) = -jti, and x = t Q , whence 


/ = .p ( x ) 


t - X 


(8-39) 


is the integral equation to he solved. From Eq. ( 8 - 38 ) the solution is 

f(x) = T dt + $ hU) (8-4o) 

it L X(t ) (t-x) 


To find X(x) in Eq. (8-23) and <$h( x ) in Eq. (8-20), get G(t) = -1 from 

Eq. (8-32), log G(b^) = rti in Eq. (8-15) so that % = - | at x = -1, 

a 2 = ^ at x = 1 with L the line -1 ^ x % 1. Thus in Eq. (8-l6), 

Ai = 0 or 1, A 2 = 0 or -1. In Eq. (8-20), with G = -1, 


whence 


P(x) = | log ( ) 

2 \ 1 + x / 

X(x) = (l+x) Al (l-x) A2 yi-^- 


Xl(x) = — I 


n/ 1 - x 2 


Xz(x) ■ VH 

X 3 (x) = 


- x 

X 


_x 

X 


X 4 (x) = V 1 - x 2 


Ai — 0 , — -1 


Ai = 0 , A2 = 0 


Ai = 1 , As = -1 


Ai = 1 j A 2 — 0 


(8-41) 

(8-42) 

(8-43) 

(8-44) 


For these four possible forms of X(x), the homogeneous solutions 0^(x) 
have the value 0 except for $^(x) = Co/'/i - x 2 (pp. 240-242, Ref. 6). 
For the case of X 4 with both ends bounded there is a restriction on F(t) 


4l 


(8-45) 


r dt = o 

_i n/ 1 - t 2 


The proper form of X(x) to he used in Eq. (8-40) depends upon the 
desired behavior of the solution at the end points. For example, in thin 
airfoil theory of aerodynamics where f(t) in Eq. (8-39) represents the 
pressure distribution, it is desired that f(l) = 0 at x = 1 to satisfy 
the Kutta condition while f(-l) may be infinite. Thus the form X 2 (x) is 
used in Eq. (8-4o) for this problem. 


Solution of the Foucault Test Integral Equation 

The Eq. (2-11) for the solid mirror is the same as Eq. (8-29) with 


A(x) = log ^ ~ X , B(x) = -jti , x = t Q , 
1 + x 


(8-46) 


and the solution from Eq. (8-38) is 


f(x) = 


F(x) log 


D 2 ( 


(l + x) + X(x) f 1 

x) D(x)JL 1 


<E> h (x) 

■ + — 

X(t)D(t)(t-x) D(x) 


F( t)dt 


D(x) = >/jt 2 + log 2 / 

V 1 + x/ 


(8-47) 


To find X(x) in Eq. (8-23) and. 0^(x) in Eq. (8-20), start with 
(see Eq. 8-32) 


G(t) = 


log 


f — 

\1 + t / 


+ jti 




= e 2i0(t) 


(8-48) 


where 


0(t) = arctan 


log 


( 


jt 


1 - t 
1 + t 



(8-49) 
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New 0(t) = 0 at t = -1, and 0(t) = jt at t = 1 so that in Eq. (8-15) 


log G(-l) = 0 , log G(l) = 2*i 


whence a x = 0 at x = -1, a 2 = 1 at x - 1. From Eq. (8-l6), Ai 
= -1 and from Eq. (8-17 ) 


T(z) = (l-z)" 1 


In Eq. (8-12) 

p(x) = -L. f 1 HZM dt 

2iti J ± t - x 

dt 

t - X 


* J -: 1 

~ [0 (t ) log (t-x)]i - 2^ 


log (t-x)dt 


(i- t2 )p + io^(Hr)' 


= log (l— x) + Q(x) 


where 


Q(x) = 


- 2 f — r 2- 

(i-t 2 ) L 2 + 


(t-x) 


To evaluate q(x) , differentiate and change variables: 

.af 

*^- 1 (l-t 2 )!* 2 + iog 2 fr 


dt 


dQ(x) 

dx 


(l-t 2 )[> + lcg i! (j^)](t-x) 


( 

■-/ 


du 


'- 00 (rt 2 +U 2 ) ' M ' e 


Or 


u 


u 


- X 
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0, 

(8-50) 


(8-51) 


where 


Take 


whence 


e n _ 1 ~ t t = -tanh - 
1 + t 2 


dt = - — sech 2 H du 
2 2 


1 - t 2 = sech 2 - 


e v = 1 ~ x 
1 + x 


dQ = 1 f° 

dx 1 + xJ_ 0 


(l+e u ) 


( e u_ e v)( rt 2 + u 2 ) 


= — f - 

1 + x -ico rt 


du 2 


du 


f 


du 


1 + x 


2 +u 2 (l+x) 2 -Loo (e u -e v )(it 2 +u 2 ) 

+ S(x) 


“I 0 

— axctan — 
n * 

-C 


= + S(x) 

1 + X 

The integral S(x) has poles in the upper plane at 


so that 


u = irt ? u = v + 2nlx 5 0 < n < °° 


S(x) = 2jri 2 Residues 

(l+x) 2 

_ 4tti 

“ (l+x) 2 

-1 


oo 

l + _l 

/, . T7- *— * 


-2rti(l+e v ) e v n=0 rt 2 + (v+2n: 

4rti £ 1 


1 " 

/■+2nirt ) 2 


l + x 1 - x 2 n _Q it 2 + (v+2nirt ) : 


(8-52) 


(8-53) 
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where Eq. (8-53) has been used to simplify the results. Use partial 
fractions and evaluate the sum to get 


S(x) = - — + ^ — [" 1 

1 + x 1 - x 2 n= o 2jr Lv + (2n+l)«i v + (2n-l)itiJ 


i— + — 2 / — 1 x 

1 + x 1 - X 2 \ v - rti/ 


_1 + 2 v + rti 


1 + x 1 - x 2 v 2 + rt 2 


Thus 


or by Eq. (8-53) 


Q 


dQ _ 2 v + rti 

dx 1 - x 2 v 2 + rt 2 


= IVHf av + *i / 

+ v 2 J rt 2 + ■ 


= log n/^+v 2 + iti 0 + C x 


and from Eq. (8-51) 


P(x) = log (l-x) + log'y^t 2 + log 2 ( ) 


+ Cl 


(8-54) 


where the rti© term has been dropped from the real P(x) ; this term 
represents the G(t) term already included in Eq. (8-24). Now from 
Eqs. (8-23), (8-50), (8-54) 


X(x) 



= C 2 D(x) 


(8-55) 


where D(x) is defined in Eq. (8-47). In this case, the only homogeneous 
solution for f(x) is a constant (see paragraph (e) in Sec. Ill) so that 
in Eq. (8-20 ) 
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3> h (x) = C^X(x) 
= C4D(x) 


( 8 - 56 ) 


The solution (8-47) now becomes 


f(x) = 


F(x) log 


\1 + X/ 


D 2 (x) 


J-i D- 


F(t) 
(t)(t-x) 


+ C 4 


D 2 (x) = jt 2 + log 2 


V 1 + x/ 


(8-57) 


which is the inversion integral solution of Eq. (2-11) . 

There may be a question about this solution (8-57) since A(x) in 
Eq. (8-46) for this problem does not satisfy the H condition at the end 
points x = ± 1, as required in Eq. (8-29). Also, this A(x) term in 
D 2 (x) in Eq. (8-57) forces the first term to be 0 at x = ± 1, which is 
a undesired restriction on f(x). To avoid these difficulties, 
log (l-x)/(l+x) can be taken as some finite value at x = ± 1, or f(x) 
can be calculated up to ± 1 + Ax. Since F(t) is known only at specified 
points, it is necessary to use numerical integration in the calculation 
of f(x) in Eq. (8-57) so that the end points can be omitted or approxi- 
mated. Since the log (l-x)/(l+x) term is almost bounded and has a rapid 
change to infinity for a very small Ax at the end points any approximation 
for the end points will have little effect on the values of f(x) at other 
points . 


Numerical Integration for Inversion Integral 

For the calculation of f(x) in Eq. (8-57) by numerical integration, 
the same procedure, except at end points., as for the collocation method 
in Sec. VI was used. The trapezoid rule was used and all log (l-x)/(l+x) 
terms in Eq. (8-57) were approximated by 

log 1 " - x . * ~ - A lt (8-58) 

1 + Xi 

as given in Eq. (6-12 ) . Thus 

Di 2 = rt 2 + (A^) 2 (8-59) 
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and 


M 

f(xi) = £ H 3 

3=1 

In matrix form, Eq. (8-60) is 

{f} = [G 0 ]{F] + C 4 [l] (8-61) 

A12 = H 2 , Ax 3 = ^ H3 , A MjM _x = - H M _x , AmjM- 2 = “ 2 Hm-2 ’ 

Goij = ; G oll = - , [1] = column of ones \ 

Dj Bf 

Df = 30.0 = D 2 , 

where A* 3 and An are defined in Eq. (6-12) except as noted above. To 
evaluate the constant C 4 , the condition f(0) = 0 was used so that 


F(x.i ) + F(xi ) 

-2d (ak) + c 4 (8-60) 

Xj - Xi 


c 4 = -■ [G 0 ,2x , j ] {Fj ) 


and 

(f) = [G]{F] 

with 

Gij = G oi j - G 0 , 2x > j > 

Oli = 0.2032 = -G mm , 
G 21 = -0.02 50 = -G m -i,m 


(8-63) 

(8-64) 


(8-65) 


The special values of if, Gxx, and G21 were arrived at by a special 
integration over the interval -0.975 to -1.0. At x = -1, the integrand 
in the integral in Eq. (8-57) becomes infinite as t -» -1. The value for 
G11, which represents this area from -0.975 to -1.0 with F = constant 
was obtained by a numerical integration from -0.975 to -0.999 and by 


J 

-1 


■0.999 


dt 


-1 -0.999 


(1+t ) log 2 ( ) log (^p) J 


= 0.132 


-1 
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where the « 2 term in ( 8 - 57 ) was dropped as small compared to the log 2 
term. The same result holds at x = 1, t = 1. For all x t points ^ ± 1, 
as t ->± 1, the integrand is 0. The integration over -0.975 to -1.00 
for these points was approximated by the Df in Eg. (8-62) and by G 21 in 
Eg. ( 8 - 65 ). 


It should be noted that the solution in Eg. (8-57) is not necessarily 
unigue. A more restricted solution, making f (x) zero at both ends, can 
be obtained by taking X(x) = constant in Eg. (8-55)* This gives the 
solution 


f(x) F(x) log (l + x) , 1 f 1 F(t)dt 

D 2 (x) D(t)(t-x) 


( 8 - 66 ) 


which differs from that in Eg. (8-57). This constant value for X(x) can 
be obtained by considering G(t) in Eg. (8-48) as a real variable expres- 
sion so that, with log(-l) = jti, 


G(t) 


w(rfr) ! 


(8-67) 


Thus, 0 = 0 in Eg. (8-49) and X(x) is constant. 

If f(x) is actually zero at the ends, then both solutions (8-57) 
and ( 8 - 66 ) may be the same. This was found to be true in the example 
of the hump discussed in Section IX. 


IX - COMPARISON OF THE VARIOUS SOLUTIONS 
FOR SELECTED EXAMPLES 


In Sections IV-VTII various solutions of the integral eguations 
(2-11 ) or (2-12) have been obtained. All of the solutions have been 
expressed in the form 


[f] = [G][F] 


(9-1) 


where [G] is a matrix that operates on a set of known values of the F(x) 
function to produce a set of values of the unknown f (x) function. The 
elements of the [G] matrix depend upon the method of solution. Since 
all of the solutions are approximate because of numerical integration or 
because of a finite functional representation or because of both of these 
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approximations, there is no reason to expect the [G] matrix to he unique 
or to expect any two solutions to have exactly the same elements in the 
[G] matrix. However, if the solutions given by the various methods are 
to agree for a broad class of F(x) functions, then there should not be 
much variation among the elements of the [G] matrices for the same 
boundary conditions on f(x). 

In order to evaluate the various solutions, several examples have 
been selected for which the exact solutions are known. The solutions 
given by the various methods are then compared to each other and to the 
exact solution for these examples. The methods of solution that do all 
the examples the best and which appear to cover the largest class of 
F(x) functions may then be considered to be the better methods of solutioi 
Also, a comparison of two selected rows of the [G] matrices is made for 
the methods in order to establish criteria as to which method of solution 
may be better. 

In order to identify all of the methods of solution, a summary 
(Table I) is presented below, which gives the form of the [G] matrix for 
each method and gives the reference equations for each matrix in [G]. 
Reference will be made to this table in the following discussion. 

To compare results given by the [G] matrices in Table I when 
used in Eq. (9-l), all calculations were made for 40 equal intervals on 
the line across the mirror, giving 4l points. Thus the [G] matrices are 
all 4l by 4l (except method 8, iteration, where the end points were 
omitted, and method 6 which was not calculated as it gives results at 
unequal Ax intervals). Examples using method 6 are given in Ref. 2. 

Example of Hump . Assume that on a line across the mirror, an 
imperfection in the form of a local hump occurs with the shape given by 


f(x) = (1 - z 2 ) 2 , z 2 < 1 , 

= 0 , z 2 > 1 

X ™ 3» 

z - , a and b are selected constants 

b 

This particular type of imperfection was used by Linfoot in Ref. 1, 
except he assumed a symmetrical case with two humps. Substitute this 
f (x) into Eq. (2-11) and integrate to get 
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Table I - Summary of [G] Matrices for All Solutions 


Method 

[G] Matrix 

Reference Eqs. for 
Matrices in [G] 

No. 

Title 

Eq. No. Formula 

1 

Power Series B 

(4-21) | [YirK] _1 [X][H] 

4-19, 4-16, 4-llb , 3-1 

2 

Power Series C 

Erj^ 

pq 

i — i 

h 

i — i 

PQ 

^pq 

i — i 
t — i 

r- 

OJ 

4-19, 4-23 

3 

Fourier Series A, 
trapezoid rule 

(5-20) 4 [S][I + K] _1 [C][H) 

nr 

5 - 18 , 5-12, 5-13, 
5-14, 5-16, 3-1 

4 

Fourier Series A, 
special rule 

(5-24) 4 [S][I + kHd] 

nr 

5 - 18 , 5-12, 5-13, 
5-14, 5-23 

5 

Fourier Series B, 
special rule 

(5-38) 4 [y][q] -1 [r][d][im] 

4- 19, 5-33, 4-8, 5-31 

5- 23, 5-35, 5-38 

6 ! 

Fourier Series B, 
Katzoff procedure 

(5-46) •- [T][Q] _1 [d] 

it 

5-47, 5-33, 4-8, 5-48 

7 

Collocation 

(6-i8) - [b rw 

6-15, 6-12 

8 i 

Iteration 

(7-16) 4 [x D^l + J + J 2 + 

* + 

7-10, 7-14, 7-11 

9 j 

Inversion Integral 

(8-65) 0. £J - G. a>t 

8 - 62 , 8-59, 6-12 



F(x) = -2z 3 + — z - (1 - z 2 ) 2 log 1-2 1 + x 
3 1 + Z 1 - X 


z 2 < 1 


= ± 2 j z = ± 1 , respectively. 


- (z 2 -l) 2 

3z 3z 3 


— + — + iog | z ~ 1 | 
z 3z 3 |z + ll_ 


> (9-3) 


z 2 > 1 


The form for z 2 > 1 has been changed from that given by the integration 
in order to obtain a form which is more accurate for large values of z. 
Graphs of Eqs. (9-2) for f(x) and (9-3) for F(x) with a = 0.6 and 
b = 0.2 are shown in Fig. 9-1. 


Because of the numerical integration in the various methods, as 
has been pointed in paragraph (b) of Section III and in Eqs. (6-4) and 
(6-5), the derivatives of the function f(x) in Eq. (9-2) should be 
examined. Now from Eq. (9-2) 


— = - — (l-z 2 ) , Z 2 £ 1 

dx b 

= 0 , z 2 > 1 , 


(9-4) 


whence f'(x) is continuous for -1 ^ x ^ 1. Thus, from Eq. (6-4), the 
integrand in Eq. (2-12) is continuous. However, the second derivative 
is discontinuous at z = ± 1, 


d 2 f _ _4_ + 12z 2 

dx 2 b 2 b 2 


z 2 < 1 


= 0 , z 2 > 1 , 


(9-5) 


whence from Eq. (6-5), the first derivative of the integrand in 

Eq. (2-12) is discontinuous at z = ±1, and the higher derivatives are 

infinite at z = ± 1. 

In order to examine the derivatives of the integrand 


g(t,x) s 

t - X 


(9-6) 
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Exact f(x),Eq.t9-2) 


Exact F(x), Eq.(9-3) 

Modification to F(x), 

See Fig. 9-2 for f(x) 





in Eq.. (2-12) for the f(x) in Eq. (9-2) in the neighborhood of z = ± 1, 
take 


z = , y = 3LZL* (9-7) 

b b 

so that 

g(t,x) = h(z,y) = (1-z ) ~ (l-y 2 ) s z 2 < 1 , y 2 < 1 , (9-8) 

b(z-y) 

, z 2 <l , y 2 ;>l , (9-9) 

b(z-y) 

= . , Z 2 >1 , y 2 < 1 , (9-10) 

b(z-y) 

= 0 , z 2 > l , y 2 > 1 . (9-H) 


Take the case of Eq. (9-10 ) and differentiate with respect to t, whence 

dg = (l-y 2 ) 2 

dt b 2 (z-y) 2 

£g = _ 2(l-y 2 ) 2 
dt 2 b 3 (z-y) 3 


d a g = /_ 1 )n+x nl(l-y 2 ) 2 

dt n b n+1 (z-y) n+1 


(9-12) 



( x ) n+i n! (l+y) 2 

b n+1 (l-y) n-1 


(9-13) 


d"g ~] = n? (l-y) 2 (9-14) 

dt*_ z = -1 b n-1 (l+y) n “ 1 


These results in Eqs. (9-12) - (9-l4) show that the higher order deriva- 
tives of g(t,x) are very large at and near z = ± 1 for all y 2 < 1. This 
implies that a lower order numerical integration rule may give better 
results than a higher order rule (see paragraph (b) in Section III) for 
this particular example. To check this implication, several different 
rules were used for some of the methods in the calculations for this 
example . 
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The matrix [F] in Eq. (9-l) was calculated for 4l values of x at 
intervals of Ax = 0.05 for -1 < x < 1 using the Eq. (9-3) with z = 5x -3 
(a = 0.6 and b = 0.2 in Eq. (9-2)). The [G] matrices as defined in 
Table I for the various methods were then used to give f t at the same 
4l points. The results are shown in Table II. 

All of the methods located the hump and all methods, except method 3, 
Fourier Series A with trapezoid rule, and method 8 using Eq. (7-5) in 
the iteration, gave a satisfactory representation of the hump. However, 
a study of the Table II indicates that some of the methods gave smaller 
maximum deviations than others. Method 9? the inversion integral, appears 
to give the best results with the trapezoid and Simpson's rules giving 
essentially the same results. Method 1, power series B, shows slightly 
larger deviations in the region of the hump for a 5 -point integration 
rule as compared to the trapezoid and Simpson’s rules. This indicates 
that even the large higher order derivatives of Eq. (9-12) - (9-14) do 
not have much effect on the accuracy of the integration rules in this 
case. However, in this method 1, the end point values are much improved 
by the higher order rules. This occurs because the spacing of the points 
is too large to give a fair representation of x n (n of order 20) when x 
approaches ± 1 at the ends. A higher order rule will naturally improve 
this representation near the ends. In method 5, Fourier series B, the 
equal Ax answer corresponds to using an interpolation matrix [IM] to 
convert the Fj values in the Table II to equal A 9 values needed in the 
calculations, while the equal A 0 answer corresponds to calculating F* 
at equal A 9 intervals directly from Eq.. (9-3) and deleting the [IM] matrix 
in [G]. The latter answer appears to be slightly better in this example. 

Since methods 3, Fourier series A with trapezoid rule, and 8, 
iteration, did not produce answers as good as the other methods in the 
hump example, and since these methods have other drawbacks already 
discussed in the description of the methods, Sections V and VII, no 
further consideration will be given to these two methods in the remaining 
examples and discussion. Also, only the equal Ax: case will be used for 
method 5, Fourier series B, in the following examples. 

Example of Modified Hump . In order to examine the behavior of the 
solutions when a large slope change occurs in F(x), a modification to 
the above hump example was made by taking F(x) = -2.0 at the points 
x = 0.40, 0.45, 0.50, 0.55 srnti 0.60 (see Fig. 9-1) • Such a modification 
not only gives an extreme change in the F(x) slope at x = 0.40 and 
x = 0.60, but also could represent a physical situation in which the 
light intensity is too weak to show a variation for F(x) < -2.0. The 
results for f(x), the mirror surface deviation, given by the various 
methods for this modified F(x) are shown in Fig. 9-2. The graph covers 
the range of 0.2 < x < 0.8. For the remainder of the mirror line, all 
methods gave small values, f(x) < 0.05, except for power series B, which 
gave larger values at x = ±1, the end points. The exact answer for 
f(x) is not known, but if the inversion integral is accepted as being 
the closest to the correct result, then the function methods, power series 
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Table II 


. Example of Hump Case (Eq. 9-2) 


Method -* 


© 

© 

© 

© 


© 

© 

CD 

(!) 

(D 

© 

® 

<D 

(9) 



f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

f 

7 

7 


f(x) 

Power Seriea 

Pcwer Seriea 

Power Seriea 

Power Seriea 

Fourier Seriea 

Fourier Seriea 

Fourier Seriea 

Fourier Seriea 

C olio cat! cn 

Collocation 

Iteration 

Iteration 

Inversion 

Inversion 

x F(x) 

Exact 

B 

B 

B 

C 

A 

A 

B 

B 





Integral 

Integral 







(Equal Ax) 

(Equal ax) 

(apecial rule) 

(apecial rule) 



Eq. (7-4) 

Eq. (7-5) 





Trapezoid 

Simpson 

5-Point 


Trapezoid 

Special Rule 

Equal AX 

Equal Afi 

Trapezoid 

Siapaon 

5th Iteration 

5th Iteration 

Trapezoid 

Siapaon 

-1.00 - .134 

0 

- .066 

- .024 

- .017 

- .015 

0 

0 

.003 

.X3 

.X5 

.X7 

.034 

.142 

.OX 

.OX 

- .95 - .138 

0 

.001 

- .001 

- .000 

.013 

.046 

-.011 

.004 

.004 

.004 

.006 

.002 

.030 

.OX 

+.001 

- .90 - .143 

0 

.014 

.016 

.020 

- .004 

.097 

-.017 

-.XI 

- ,X2 

•X5 

.006 

.020 

.021 

-.OX 

+ .0X 

- .85 - .148 

0 

.022 

.018 

.017 

- .003 

.045 

-.005 

.005 

.X5 

•X5 

.006 

.029 

.037 

-.ox 

+.0X 

- .80 - .153 

0 

.003 

.003 

.002 

.006 

.074 

-.014 

.003 

•X3 

.005 

.X6 

.028 

.040 

-.ox 

.ox 

- .75 - .159 

0 

.008 

.006 

.011 

- .003 

.104 

-.021 

-.X2 

-.004 

•X5 

.006 

.026 

.04C 

-.ox 

.ox 

- .70 - .165 

0 

.021 

.019 

.022 

- .007 

.075 

-.012 

.OX 

-.XI 

.005 

.006 

.024 

.039 

-.coo 

.ox 

- .65 - .171 

0 

.018 

.017 

.017 

.000 

.043 

-.005 

.006 

.006 

•X5 

.006 

.022 

.038 

-.ox 

.ox 

- .60 - .178 

0 

.007 

.007 

.006 

.005 

.058 

-.010 

.X6 

.006 

.X5 

.006 

.020 

.035 

-.ox 

.ox 

- .55 - .186 

0 

.004 

.004 

.005 

.000 

.094 

-.020 

-.XI 

-.XI 

.X5 

.00 6 

.018 

.033 

-.ox 

.ox 

- .50 - .195 

0 

.012 

.on 

.014 

- .007 

.104 

-.020 

-.004 

-.X5 

.004 

.006 

.016 

.030 

-.ox 

.ox 

- .45 - .204 

0 

.021 

.U9 

.021 

- .007 

.076 

-.012 

.XI 

-.XI 

.X5 

.006 

.014 

.027 

-.ox 

.ox 

- .40 - .215 

0 

.020 

.018 

.019 

.001 

.045 

-.004 

.X7 

.X7 

.004 

.006 

.013 

.024 

-.ox 

.ox 

- .35 - .226 

0 

.010 

.010 

.009 

.006 

.044 

-.006 

.X8 

.X9 

.X5 

.006 

.011 

.021 

-.ox 

.ox 

- .30 - .239 

0 

.003 

.003 

.003 

.002 

.074 

-.016 

.002 

.002 

.004 

,005 

.009 

.018 

-.ox 

.ox 

- .25 - .253 

0 

.005 

.005 

.007 

- ,00c 

.106 

-.023 

-.005 

-.006 

.006 

.X7 

.007 

.017 

-.ox 

.ox 

- .20 - .269 

0 

.015 

.013 

.017 

- .010 

.108 

-.020 

-.004 

-.006 

.X3 

.004 

.006 

.012 

-.ox 

.ox 

- .15 - .287 

0 

.023 

.021 

.023 

- .004 

.075 

-.009 

.X3 

,X2 

.X7 

.007 

.004 

.009 

-.ox 

.ox 

- .10 - .308 

0 

.021 

.020 

.020 

.006 

.029 

-.001 

.010 

.011 

.002 

.003 

.003 

.007 

-.ox 

.ox 

- .05 - .333 

0 

.010 

.010 

.009 

.009 

.000 

-.004 

.X9 

.011 

.008 

.008 

0 

0 

-.ox 

.ox 

0 - .361 

0 

0 

0 

0 

0 

.00c 

-.017 

0 

0 

0 

0 

-.002 

-.002 

0 

0 

.05 - .396 

0 

.000 

- .001 

.002 

- .012 

.037 

-.028 

-.X9 

-.011 

.011 

.010 

-.003 

-.006 

.ox 

-.ox 

.10 - .437 

0 

.013 

.011 

.015 

- .014 

.070 

-.024 

-.X7 

-.010 

-.003 

-.004 

-.005 

-.008 

.ox 

-.ox 

.15 - .488 

0 

.029 

.027 

.030 

.000 

.087 

-.006 

,X7 

.005 

.015 

.012 

-.006 

-.010 

.ox 

-.coo 

.20 - .554 

0 

.032 

.030 

.030 

.016 

.081 

.009 

.018 

.020 

-.008 

-.010 

-.008 

-.012 

,ox 

-.ox 

.25 - .642 

0 

.012 

.011 

.008 

.013 

.061 

-.000 

.010 

.013 

.022 

.017 

-.010 

-.013 

.ox 

-.XI 

.30 - .765 

0 

- .017 

- .020 

- .021 

- .016 

.044 

-.032 

-.015 

-.015 

-.016 

-.021 

,010 

-.013 

.ox 

-.XI 

.35 - .956 

0 

- .020 

- .024 

- .019 

- .035 

.056 

-.048 

-.022 

-.030 

.039 

.031 

-.015 

-.009 

-.002 

-.005 

.40 -1.333 

0 

.055 

.053 

.062 

.023 

.131 

.019 

-.051 

.036 

-.030 

-.038 

-.003 

.027 

.X5 

.004 

.45 -2.214 

.191 

.247 

.249 

-.256 

.210 

•307 

.216 

.248 

.233 

• 193 

.2X 

.184 

.173 

.200 

.206 

.50 -2.653 

.562 

.542 

.547 

.545 

.535 

.587 

.526 

.547 

.543 

.573 

.573 

.532 

.439 

.561 

.563 

.55 -2.338 

.879 

.850 

.854 

.843 

.858 

.900 

.836 

.841 

.853 

.890 

.893 

.845 

.703 

.874 

.874 

.60 -1.386 

1.000 

1.020 

1.019 

1.010 

1.014 

I.O91 

.987 

.977 

.999 

1.013 

1.020 

.970 

.830 

.993 

.997 

.65 - .112 

.879 

.925 

.922 

.926 

.089 

1.006 

.872 

.862 

.880 

.883 

.884 

.862 

.768 

.871 

.871 

.70 1.059 

.562 

.576 

.579 

.595 

.544 

.642 

.540 

.546 

.549 

.572 

.583 

.568 

.543 

.556 

.560 

.75 1.656 

.191 

.179 

.185 

.196 

.188 

.228 

.186 

.211 

.201 

.165 

.171 

.220 

.249 

.190 

.193 

.80 1.333 

0 

- .000 

- .008 

- .016 

.009 

.056 

-.000 

.026 

.016 

.018 

.017 

.034 

.048 

.010 

.015 

.85 .956 

0 

.046 

.030 

.019 

- .005 

.122 

-.014 

-.X3 

-.X3 

.007 

.008 

.026 

.024 

-.X3 

-.X3 

.90 .765 

0 

.007 

.029 

.040 

.001 

.080 

-.010 

.004 

.006 

•X5 

.007 

.XI 

.040 

-.001 

-.ox 

.95 .642 

0 

- .019 

- .029 

- .027 

- .002 

-.019 

-.016 

.OX 

-.OX 

•X5 

.006 

-.096 

-.191 

.002 

-.X3 

1.00 .554 

0 

- .272 

- .104 

- .076 

- .003 

0 

0 

.002 

.002 

.005 

.X7 



.ox 

+ .0X 

Max. + deviation 

+ .056 

+ .058 

+ .065 

+ .027 

+ .131 

+.025 

+ .057 

+ .042 

+ .039 

+ .031 

+ .034 

+.142 

+.018 

+.015 

Max. - deviation 

- .272 

- .104 

- .076 

- .035 

-.019 

-.048 

-.051 

-.030 

-.030 

-.038 

-.096 

-.191 

-.008 

-.000 


See Fig.(9-I) For Modified F(x) 
a Power Series B (trapezoid rule) 
v Power Series C 
□ Fourier Series A (special rule) 

0 Fourier Series B (equal Ax) 

-o Collocation (Simpson's rule) 

— Inversion Integral (trapezoid rule) 


/ 


s;. 


Fig. 9-2. f(x) for Modified Hump Case 





and Fourier series, appear to oscillate about the inversion integral 
curve and thus represent it approximately. The representation appears 
to be as good as could be expected with twenty constants in the function 
solutions. The collocation solution appears to represent the inversion 
integral solution with a high frequency large amplitude oscillation. 

Since the collocation method has no smoothing function effect in it, it 
may be expected to oscillate on functions with large slope changes or 
with discontinuities. Further discussion of this deficiency of the 
collocation method in given below in the G matrix paragraph. 

Example of Straight Line . As pointed out in paragraph (f) of 
Section III, the method of solution should properly solve the case of a 
constant F(x), which produces a rotation of the mirror or a straight 
line deviation across the mirror. To examine this case, [F] in Eq. (9-1) 
was taken as Fi = -40.0 at all 4l points so that the exact slope of f(x) 
is 20.0 with f(-l) = -20.0, f(x) = 20.0, and f 21 = 0, fa 2 = 1.0, etc. 
Methods 1, 2, 4, 5j 7, 9 ti Table I were used with the following results. 

Method 1, power series B with trapezoid rule, gave a large oscilla- 
tion with deviations of order of ± 2.0 from the straight line and with 
+ 38.78 instead of + 20.0 at the ends. As pointed out above in the 
discussion of the hump example, this failure of the Power series B method 
is due to the poor representation of x n near the ends. This affects the 
G matrix, which shows large effects when f (x) has its largest value at 
the ends. See discussion of the G matrices below. 

Method 2, power series C, gave the straight line to five significant 
figures at all 4l points. 

Method 4, Fourier series A with special rule, did not give the 
straight line. It has large oscillations with deviations as large as 
± 10.0, with ±20.0 deviations at the ends. This failure of Fourier 
series A is due to the restraints imposed by the sin nx functions being 
zero at the ends . This restraint gets into the G matrix, which means 
that for f(x) odd and not zero at the ends, the method will give poor 
results. In this case, the series must represent a function with a 
discontinuity at both ends. See discussion of the G matrices below. 

Method 5, Fourier series B with special rule and equal Ax’s, gave 
the straight line to four significant figures at all 4l points. 

Method 7, collocation with Simpson's rule, gave the straight line 
to at least eight significant figures (only eight digits were printed 
out) at all 4l points. 

Method 9 5 inversion integral with trapezoid rule, gave the straight 
line to four significant figures at all 4l points. 

The above three examples of the hump, the modified hump, and the 
straight line show that all of the methods except possibly the inversion 
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integral, appear to have deficiencies. Some methods do better on one 
type of function than on other types. Althougn power series B and 
Fourier series A did satisfactorily on the hump and agreed with the other 
function methods on the modified hump, they failed completely on the 
straight line. The collocation method did very well on the hump and the 
straight line, but did poorly on the modified hump. The power series C, 
the Fourier series B, and the inversion integral methods appear to be 
the most consistent in the three examples considered. In order to 
determine which may be the better methods for all types of possible 
functions f(x) and F(x), an examination of the G matrices was made for 
the various methods. 

The G Matrices . To study the behavior of the G matrix in Eq. (9-1) 
for the various methods of solution, graphs of two typical rows of the 
G matrices were constructed. Figure 9-3 shows row 16, which gives f 16 
at x s -0.25, and Fig. 9-k shows row 33, which gives f33 at x = 0.60. 
Methods 1, 2, 4, 5? 7, 9 Table I are shown in Figs. 9-3 and 9-^* 

Method 1, power series B with trapezoid rule, method 2, power series C, 
and Method 4, Fourier series A with special rule are too close to dis- 
tinguish on the graph, except at the ends. For row 16, at x = -1, 
point 1, power series B gives 0.0077, power series C gives 0.0000, and 
Fourier series A gives 0.0538; for row 16 at x = 1, pt. 4l, power series B 
gives 0.0191, power series C gives 0.0001, and Fourier series A gives 
0.0520. 

Since these differences in the end values occur in all the rows of 
the G matrices for methods 1, 2, and 4, it appears that these large end 
values in power series B and Fourier series A are the cause of the failure 
of these two methods to give the proper results in the straight line 
example above. Since the power series C has small end values in the 
G matrix, agrees with the power series B and Fourier series A elsewhere, 
and gives the proper straight line, it would appear to be the best of 
these three function methods. 

In Figs. 9-3 and 9-^+, the Fourier series B deviates slightly from 
the other three function methods and has small values at the ends 
(Gl6»l = 0 .0003, Gx 6,41 = -0.0028). On the whole, it should give results 
similar to power series C. 

In Figs. 9-3 and- 9-^+, if appears that the function methods are 
actually representing the inversion integral curve as well as they can 
with twenty constants . The inversion integral G matrix, -which is limited 
in accuracy only by the numerical integration for each row, should be 
the best overall matrix (except possibly the first row and last row). 

The peak values in the rows of the G matrix occur on either side 
of the origin (x =0), which is the restraint point for the solution to 
make f(0) = 0, and on either side of the f* point being calculated 
(Gi6,i5 and. G 16 ,i 7 in Fig. 9-3 for row 16 giving f 16 ). Since the form 
of the integral Eq. (2-ll) gives large values in the integrand for 
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t near x, the point at which f (x) is being calculated, it is evident 
that these adjacent values should be emphasized in the solution. Of 
necessity, the function methods must smooth out these peak values and 
fail to give proper emphasis to them. For these reasons, the inversion 
integral method is considered to be the best method for all types of 
possible functions. 

The entire G matrix for the inversion integral method using the 
trapezoid rule for the numerical integration is given in Appendix C. 

It may be noted that the inversion integral G matrix with Simpson's rule 
has slightly larger peak values than those shown in Figs. 9-3 and 9-4 
for the trapezoid rule. 

The collocation method shows a high frequency oscillation, with 
large amplitudes, about the inversion integral curve in Figs. 9-3 and 
9-4. It is evident that if F(x) has a sharp change or is discontinuous, 
then the large values in the collocation G matrix can magnify the dis- 
continuity and give incorrect values for f(x). Thus, the collocation 
method, in spite of its performance on the straight line and the hump 
example, it is not a good method for all types of functions, particularly 
those with discontinuities. 

Evaluation of the Various Methods . On the basis of the above 
examples and examination of the G matrix, an evaluation of the nine 
methods of solution listed in Table I can be made. Although all the 
methods can solve certain problems, the best methods are those that can 
solve the largest class of problems. Since F(x) must be measured at 
points, and f(x) may represent any conceivable deformation of the mirror, 
it is essential that the method of solution handle almost any type of 
function, whether continuous or not. From this viewpoint, the following 
evaluation is made: 

1. Method 9s "the inversion integral, is the best overall 
method. It will handle most any case. 

2. Power series C, method 2, and Fourier series B, 
methods 5 and 6, are the best function methods. They 
will handle most any case. 

3. Collocation, method 7s is simple and satisfactory for 
smooth functions, but is no good on discontinuous 
functions. A violent oscillation in the solution 
indicates the method has failed. 

4. Power series B, method 1, and Fourier series A, 
methods 3 and- 4, are unsatisfactory because of end 
points representations and restrictions. 

5o Iteration, method 8, is unsatisfactory because of end 
point restraints and numerical integration difficulties. 

Up to this point the discussion and evaluation of the various 
methods has been concerned with the solid mirror. Section X takes up 
the Cassegrain mirror, or the mirror with a central hole. 
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X - SOLUTIONS FOR CASE OF A CENTRAL HOLE 


The integral equations to be solved for the Cassegrain mirror are 
given in Eqs. (2-13) or (2-15). Take Eq. (2-15) in the following form 


F(x) 


c *n 


- f ( x ) at 

t - x 


i 


( 10 - 1 ) 


where R is defined in Fig. 2-3. Theoretically, this Eq. (10-1) can be 
solved by any of the methods described in Sections IV-VIII. However, 
on the premise that the evaluation of the methods for the solid mirror 
in Section IX holds for the Cassegrain mirror, only the better methods 
will be used for this case. The power series C, collocation and inversion 
integral methods of solution, as modified for the hole case, are described 
below. The Fourier series B, method 6 in Table I, is described in Ref. 2 
for the hole case, together with examples. It will not be covered here. 

Power Series C for Hole. As in Section IV, take 


f(x) = £ A J* J 

0=1 


in Eq. (lO-l) and integrate to get 


( 10 - 2 ) 


whence 


J J 

j=l k=l 

[1 - Rj- k+1 ] [1 - (-l)J- k+1 ]x k_1 

(10-3) 

0 - k + 1 

F( Xi ) = Ft = 

J J 

£ A jBj(xi) = £ B U A J 

0=1 0=1 

(10-4) 

_ . f [1 - 

E i-k+l] [ x _ (_i)J-k+i] Xi k -l 

(10-5) 


0 - k + 1 


These Eqs. (10-4) and (10-5) correspond to Eqs. (4-22) and (4-23) so 
that from Eqs. (4-26) and (4-28), 

[f] = [G][F] ( 10 _ 6 ) 

[G] = [Y][b!b]“ 1 [B]1 (10-7) 

where the elements of [Y] and [B] are given in Eqs. (4-19) and (10-5), 
respectively. 
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V 


Results using power series C for several examples with, a hole are 
given below and compared to the results of the other methods. 

Collocation Method for Hole . For 2R as the width of the hole on 
the particular line on the mirror (see Fig. 2-3), take (m/ 2)-1 = 20 
equal intervals on each side. This gives 21 points on each side, or 
M = 42 total number of points, including four edge points. Take 


whence 


P = Ax = 


g(l-R) 
M - 2 


M points. 


xj = -1 + P(j-l) > 0=1, | 

xj = R + p(j - 1 - M \ o = - + 1 , 

' 2 ' 2 

As in Section VI, take Eq. (lO-l) in matrix form 

[A][f] = - [F] 

where, by using Eqs. (10-9) and (6-12), 


( 10 - 8 ) 


(10-9) 

M 1 


( 10 - 10 ) 


A i,t-1 = “ Hi-1 - ^ , A ljl+1 = Hi+2. + ^ Hi , 


i^l , — , — +1 , M , 

2 2 

Al 2 = 2Hi + H 2 , A 13 = - i (Hi - H 3 ) , 

A (i> I - 2 ) = I [ H (|) - (|- 2 )] • 

A (f + 1 ■ l +3 ) = - IKI + x ) - K! + 3 )> 

A (I> I- 1 )- -KfMf- 1 ) > 

A (| + 1 , | + 2) = 2h(| + l) + H (| + 2) , 

A(M,M-2) = i [H(M) - H(M-2) ] , 

2 

A(M,M-1) = - 2H(M) - H(M-l) , 


H.,P 


Aij = -141 


xj - Xi 


M 


, except above and i / j 


A ii = " H A lJ 
0=l»0^i 


\ (10-11) 
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Now, from paragraph (e) of Section III, there are two homogeneous 
solutions for the hole case. Thus two points must be fixed on the mirror 
in order to obtain a solution by the collocation method. Any solution 
thus obtained will include the homogeneous solutions (3-6) and (3-7) or 


f h = c i + — (10-12) 

x 

The value of C x is of no consequence as it represents a translation of 
the reference for f(x), but the value of C2 should be as small as possible 
to keep the apparent distortion C2 /x small. In most cases the smallest 
value of C 2 will be given by fixing the two points at the edges of the 
hole. However, if f(x) is larger at the edges of the hole than at the 
ends, then the ends would give a smaller C 2 . Note that C2 = 0 for f(x) 
even. 


Use the conditions 

f(|) 3 f(-R) = 0 

f(| + 1^ = f(R) = 0 

so that, corresponding to Eq. (6-15 ), 

[B] = [A] t [A] + [S] 

S*j = 0 5 except 

s (| M), s (! + 1 M +1 ) =1>0 I 

This gives (see Section VI ) 

[f] = [G][F] 

[G] = - [B] _1 [A]^ 


(10-13) 


(10-14) 


(10-15) 


where [A] is defined in Eq. (10-11) and [B ] in Eq. (10-14). 
Results using the collocation method are given below. 
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Inversion Integral for Hole Cage . The integral Eq. (2-13) for the 
Cassegrain mirror is the same as Eq. (8-29) with 


A(x) = log B (x) = - rti 

ll + x x - Rl 

x = t D , L = lines -1 to -R and R to 1 

and the solution is given by Eq. (8-38) as 

f (:0 _ M log 2_t*| + lisl j F(t)dt 

D 2 (x) ll + x x - Rl D(x) l X(t)D(t)(t- 

. ®hM 


(10-16) 


D(x) 

D 2 (x) = it 2 + log 2 |i— ^ 

ll + 


(10-17) 


,2 |1 - _X X + R 
x x - R 


To find X(x), Eq. (8-23), and $ n (x), Eq. (8-20 ), to use in the above 
Eq. ( 10-17 ), start with Eq. (8-32) for G(t), whence 


where 


G(t ) 


log 


log 


1 - t t + R 
1 + t t - R 
1 - t t + R i 

1 + t t - rI 


+ iti 


rti 


e 2i0(t) 


0(t) = arctan 



jt 

1 - t t + R l 
1 + t t - R - 


(10-18) 


(10-19) 


From Eqs . (8-15), (8-16), (10-18), and (10-19), values of 0 k , 
log G(b k ), a k , and A k , at the edge points can be obtained as follows: 


x 0 log G a k 7\ k 

-1 0 0 0 0 


-R it 2«i 


1 -1 s 


R 0 0 0 0 

1 jt 2iti 1 -1 


( 10 - 20 ) 
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Thus, Eg.. (8-17) gives 


T(z) = (z+R)“ 1 (1-z)“ 1 (10-21) 

From Eqs . (8-12) and (10-18) 

P(x) = - / 5 (10-22) 

rt L t ' x 


where L is -1 to -R and R to 1 and 0(t) is given by Eq. (10-19). By 
analogy to the integration of the corresponding Eq. (8-51) for the solid 
mirror with the result in Eq. (8-54), the solution for Eq. (10-22) was 
deduced as 


P(x) = log 



+ log 2 


1 - x x + 
1 + x x - 



+ C x 


(10-23) 


This result was verified very closely by a numerical integration of 
Eq. (10-22) using the trapezoid rule and the [A] matrix in Eq. (lO-ll). 
That is. 


[P] = - — [A][0] , 42 points , (10-24) 

jt 


where An =0. 

Put Eqs. (10-23) and (10-21) into Eq. (8-23) to get 
X(x) = C 2 (l-x)“ 1 (x+R)“ 1 (l-x) (x+R) x - 1 D(x) 

= c 2 , (10-25) 

x 


where D(x) is defined in Eq. (10-17). In this case, the homogeneous 
solution for f(x) is given by Eqs. (10-17) and (8-20) as 


fh( x ) 


$h( x ) 

D(x) 


DU) 

xD (x) 


(C3X 


+ C4) 


= C 3 


+ C4 


x 


( 10 - 26 ) 
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which agrees with the homogeneous solution given in paragraph (e) of 
Section III. Equation (10-17) now becomes 


f(x) 


IM- i og + i / + c 3 

D 2 (x) 1 + xx-R xl D 2 (t)(t-x) 


+ 


C4 

x 


(10-27) 


D 2 (x) = « 2 + log 2 


1 - x x + R 
1+xx-R 


( 10 - 28 ) 


which is the general inversion integral solution of Eq. (2-13) • 

It should be noted that this solution in Eq. (10-27) is not unique . 
Other forms for X(x) in Eq. (10-25 ) are 


X(x) = CaD(x) and X(x) = C 2 


( 10 - 29 ) 


■whence other solutions are 


f(x) = Jgl log 13^* SJL£| + / + c 3 + Sa , (10 

D 2 (x) ll + X X - Rl r D 2 (t)(t-x) X 


■ 30 ) 


f(x) = log I 1 - - - x — - R | + - / f^ dt + C T + ^ (10-31) 

D 2 (x) ll + x x - Rl B(x) t, D(t)(t-x) x 


D(x) i D(t)(t-x) 


It was found that the particular integrals in the solutions (10-27), 
( 10 - 30 ), and ( 10 - 31 ) behaved as if two restraints were present, or 
f(±a) = 0, for x = ± a points on the mirror. The value of a is different 
in each solution for the same function. For these reasons the mirror 
was restrained at the edge points, x = ± R, or f(+R) = 0, which corre- 
sponds to the restraints used in the collocation method. Note that by- 
replacing t by t - x + x in the integrand in Eq. (10-27), the solution 
in Eq. (10-30) is obtained with a term (C 9 /x) added. Since the solution 
in Eq. ( 10 - 31 ) makes f(l) = f(-l) = 0, it is more restricted. However, 
for the ’'hump" example given below, all three solutions are essentially 
the same. 

The numerical integration in Eq. (10-27) was carried out in the 
same manner as for the solid mirror, Eqs. ( 8 - 58 ) - ( 8 - 65 ), where all log 
terms are approximated by -An in Eq. (10-11). Thus 

f( Xl ) = y Ma*). r xjF&i) + *%11 + c 3 + s* ( 10 - 32 ) 
pi - Xi L XjDJ Df J x t 
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Df = jt 2 + (An) 2 


[f] = [Go][F] + C 3 [l] + Cj— 

Lxi. 


G °i<3 = J G oii = “ ^2" 

XiDf Df 


[1] = column of l's , I — I = column for — 

LXiJ Xi 


(10-33) 

(10-3.4) 


(10-35) 


where from Eq. (10-11) and from special integrations for the edge values 
(see Eqs. (8-62) and ( 8 - 65 ) and discussion after Eq. ( 8 - 65 )): 


Aij = 


JkZ 

Xj " X! 


i / <i > 


Ai , 1-1 = - Hi-! - | Hi , i ^ 1, ^ + 1, 

^ 2 

Ai,i«. = Hi +1 + | Hi , i 4 | > M, 

M 

Aii = - S, Ai 3 > i ¥ 1, f, f + 1, M 

j=l,j^i 2 2 


Df = D 2 = 28.0 , D 2 (|) = d£? (! + l) = H.O , 

0'cfl.i = “ GomM = 0.2060 , 

G °(f,|)= " Go (l + ls f + 1 ) = -°* 3180 ’ 

G d (2,l) = - G 0 (M-1,M) = -0.0283 , 

G 0 (| - 1,|^ = -Go(|+ 2,| + 1^ = 0.0432, 

G °(l " X, 1 + X ) = " G °(f + 2, l) = -°- 0031 


>(10-36) 


J 


> 


(10-37) 
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As noted above in the collocation method, the smallest value of 
C 4 /x in most cases will occur if the two points at the edges of the hole 
are fixed. Thus from Eqs. (10-13) and (10-34) 


0 = [ G o(|,d)[ F jJ + °3 - TC 

0 " [ G =(l +1 ' 3 )] [ F j] + ° 3 + f 

- 1 N!’ 3 ) + g °(! +i ^)]N 
|[ G °(l 3 )' G °(f +1 ’ 3 )]W 

Now, Eq. (10-34) becomes 

[f] = [G][F] 


whence 

C 3 = 
C 4 = 


(10-38) 


(10-39) 


a ‘< - q o» -K 1 - ^) g °(I»j) -l( i + ^) a »(l +i > 3 ) (i °- 4o> 


where G 0 ^ is defined in Eqs. (10-35) - (10-37). 

Results for the inversion integral method are given below and 
compared to the other methods. 


Compariso n of Solutions for Cassegrain Mirror for Selected Examples . 
The three methods, power series C, collocation, and inversion integral, 
described above were applied to several examples and the results compared. 


1. The example of a hump used above in Section IX for the solid 
mirror was also used for the Cassegrain mirror. In this case 


(1-z 2 ) 2 

, z 2 <; 1 


0 

, z 2 > 1 


x - a 
b 

, a = O.65 , b = 0.175 , 

(10-41) 

0.125 
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whence 


1-zl+xx-R 


F(x) = - 2z 3 + ^ z - (l-z 2 ) 2 log 


, z 2 < 1 


l + zi_xx + R 


=± 5 • , " tl ■ 


(10-42) 


■ S ■ + 3? - <z2 - 1) 


2 T ^2 


_ 2 _ 

3z £ 


+ -Z-S + l0 § 


(rfr) 


, z 2 > 1 


This F(x) was calculated at 42 points (20 equal intervals on each side) 
(R = 0.125) and the G matrices used to calculate f(x) at the same 42 
points. The results were similar to those in Table II for the solid 
mirror with the following maximum deviations from the exact answer in 
Eq. (10-41): 


P.S.(C) Collocation 


Deviation +O.O 65 
-0.048 


+0.033 


- 0.030 


Inversion 

Integral 

+0.017 

-0.009 


(10-43) 


It should be noted that the restraints in the power series C solution 
are that f(0) = 0, which is not on the mirror, and f(x) be continuous 
across the hole. This means that this solution may have to be translated. 
In fact, a translation of + 0.08 was used for the power series C solution 
to get the best results at the hump. Even so, the power series C solution 
deviated more than the other two solutions. Since f(±R) = 0 was true 
in this example, it appears that the restraints of f(±R) = 0 used in the 
collocation and inversion integral solutions had essentially no effect. 

2, However, in the example of F(x) = -C = constant, the restraints 
will affect the f (x) solution. With no restraints or no homogeneous 
solutions 


F = - C gives f(x) = 


2(1-R) 


(10-44) 


On the other hand, with restraints f(+R) = 0, and 


* U) = TTT-r X + c x + £2 
2 (l-R) x 


(10-45) 
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it follows that 


whence 


Ci = 0 , C 2 = - 


R" 


2(1-R) 


2 (l-R) \ x / 


(10-46) 


(10-47) 


The three methods were used for F(x) = -35 or C = 35, and R = 0.125, 
whence 


f = 20x in Eq. (10-44) , (10-48) 

f = 20x — in Eq. (10-47) (10-49) 


The power series C method gave f(x) in Eq. (10-48) to three places. 

The collocation method gave f(x) in Eq. (10-49) to four places. The 
inversion integral method gave Eq. (10-49) to three places. 

The above example of F(x) = constant has other implications. In 
paragraph (f) of Section III, it was pointed out that the reference scale 
for F(x) may be arbitrary so that a constant change in F(x) may be only 
an apparent rigid body rotation of the mirror and not a real distortion 
on the line. If an apparent rigid body rotation is occurring, then the 
collocation and inversion integral methods will show not only the rotation 
but also an apparent distortion of the form C 2 /x, due to the two fixed 
points on the mirror surface used in these methods . The power series C 
and Fourier series B methods will show only the rotation, as these methods 
do not restrain points on the mirror. If the function methods were set 
up to give two restraints, they would also show the apparent C 2 /x term 
for a constant change in F(x), whether real or apparent. 

On the other hand, if a real C 2 /x type of distortion occurred in 
the Cassegrain mirror on some line, it is not evident that any of the 
methods of solution can detect it, since F(x) = 0 in this case. Appar- 
ently, some other method of measurement would have to be used to detect 
such a distortion. However, it seems unlikely that a distortion of the 
form C 2 /x would occur on a mirror with R relatively small. The inherent 
large stiffness in the circumferential direction would tend to restrain 
the mirror from such a distortion. 

3. The three methods were applied to the example of F(x) = - 35.0 
for x S -0.125 and F(x) = 35.0 for x § 0.125, a, case similar to one used 
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by Katzoff in Ref. 2. In this case, f(x) is an even function and the 
C 2 /x term in the collocation and inversion methods is zero. After 
translating the power series C solution, it was found that all three 
solutions agreed within ± 2 in the third place. The shape of the curve 
is the same as that in Ref. 2. 

4. The three methods were applied to the example of 

F(x) = ^ , f(x) = \ (10-50) 

X X 


Here f(±l) = 1.0 and f( ±0.125) = 64.0. After translation of all three 
solutions to give f(±l) = 1.0, the maximum deviations from 1/x 2 for the 
42 points were +0.1 and -0.4 for inversion integral, +0.2 and -2.1 for 
collocation, +4.9 and -1.2, with considerable oscillation, for power 
series C . As might be expected the power series C did poorly on this 
negative power function. On the other hand, since f(x) is even in this 
case, power series C was able to translate and do a reasonable represen- 
tation of f(x). 

5. To further investigate the representation of a negative power 
deviation by power series C, the example 

F(x) = ^ , f(x) = (10-51) 


was solved. Since this is an odd function, and power series C must be 
0 at x = 0, it is a difficult function to represent, going from 
f(l) = 1.0 to f(o.l25) = 512.0. As expected, it did a very poor job, 
giving f(l) = 126, f (0.125) = 472, f (0.16875) = 510 instead of 208.1, 
f(0.3) = 370.0 instead of 37.0, etc. Thus, it would appear that the 
function methods, both power series and Fourier series, cannot handle 
odd negative power type deviations. Of course, if desired, negative 
power terms could be added to the power series methods so that they could 
represent these for the Cassegrain mirror. 

This example 5 in Eq. (10-51 ) takes the form 

F(x) = , f(x) = i - — (10-52) 

X X ' 3 X 


when the points x = +R are restrained in the collocation and inversion 
integral methods. In this particular example, the C 2 = -64 would be 
much smaller (C 2 = -l) if the end points x = ± 1 were restrained instead 


72 



of the points at the edges of the hole. However, since the G matrix for 
the inversion integral was set up for the edges x = ± R restrained, 

Eq. (10-52) was compared to the inversion integral solution. The maximum 
deviation was at x = ±0.16875, where f(x) = +183 instead of +172 in 
Eq. (10-52). By changing C 2 = -64 to C 2 = -65-9} the maximum deviation 
was less than +1.0. This indicates there are slight errors in the 

approximations for the elements an<i ®o(lf + ■*•>! + ^ 

in Eq. (10-37), which are introduced as a C/x error through Eq. (10-40) 
for the restraint conditions. 

The conclusion is that example 5 S Eqs. (10-51) and (10-52), is a 
rough shape to calculate numerically. It is also -unlikely that the 
mirror could experience such a deformation. The power series can not 
get f(x) in Eq. (10-51), hut the inversion integral gets the f(x) in 
Eq. (10-52) satisfactorily. 

G Matrices for Cassegrain Mirror . To better understand why the 
power series C solution was very poor in some of the examples considered, 
the G matrices for the three methods were examined. Figure 10-1 shows 
a graph of row 33 for the three methods. The inversion integral row 33 
in Fig. 10-1 is nearly the same as that in Fig. 9-4 for the solid mirror 
case. The collocation row 33 in Fig. 10-1 is quite different from Fig. 
Fig. 9-4, having much less oscillation. In fact, the collocation graph 
follows the inversion integral quite closely, except near column number 
33 where it overemphasizes the local effect for point 33. Apparently, 
the restraint of two points in this case as compared to one point at the 
origin in the solid mirror case produced a damping effect on the colloca- 
tion G. matrix. The power series C row 33 in Fig. 10-1 is radically 
different from Fig. 9-4, having much larger oscillations and failing to 
follow the inversion integral around column 33. This difference, which 
occurs in all the rows of the G matrix, explains why power series C does 
a poor job on certain types of functions in the hole case. 

The reason why the power series C method has the large change in 
the G matrix from the solid mirror to the Cassegrain mirror appears to 
be in the restraints for the hole case. No points are fixed on the 
mirror. Only the origin, which is off the mirror, is fixed. With no 
points being specified in the hole, the power series C functions appear 
to be insufficiently defined to produce the proper G matrix. Also, it 
appears that the frequency of the oscillation in Fig. 10-1 is determined 
by the hole size, rather than the mirror size. This may indicate that 
the number of points and the value of R, or the ratio between R and 
interval size, affect the behavior of the G matrix for the power series 
C method as well as for other function methods. 

Evaluation of Methods for Cassegrain Mirror . From the above 
examples and discussion of the G matrices, an evaluation of the three 
methods of solution used for the hole case is as follows: 
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For a specified tj(x) take 


g(xi) = at^xt) , 

Si = sin g (xi) , 

Cj = cos g(Xi) , 

[T t ] = [A] [Si ] , and 

[U t ] = [A] [C t ] , 


(11-3) 


where the matrix A is the numerical integration matrix in Eq. (10-36) 
with rows 1, M/2, (M/2) + 1, M omitted. Since the integrals become 
infinite at the edge points, F(x) at x = ±1 and x = ± R cannot be 
calculated. However, these points are included in the numerical inte- 
gration for F(xi) at all interior points. At x t , Eq. (11-2) gives 


F(xi) = Ui 2 - log 2 I 1 “ x i , x j + R 

ll + Xi Xi - R 

+ 2jtSiUi - 2jtCiTi + Tf (11-4) 

Assume g(x) to have the form of Eq. ( 9 - 2 .) as used in the linear case 
above 


g(x) = p(l-z 2 ) 2 , z 2 £ 1; 


0 

x 


z 2 > 1; and 


a = 0.6 


z 

R = 5/41 , ax = 2/4l , 


b = 0.2 


(H-5) 


where p is a factor for the magnitude of the error. Note that 


or 


p = (2rt)(max. error in wavelengths) (ll-6) 

Imax = max * error in wavelengths = p/2jt . (ll-7) 
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For the linear case used above p = l/2jt and this case was used in 
Eq. (11-4) to compare to the linear results. Figure 11-1 shows the 
comparison for F(x) using 38 -point numerical integration. The linear 
numerical integration and the nonlinear numerical are quite close. 

Figure 11-2 shows results of Eq. (ll-4) for several values of p; i.e.. 


1 

2 * ’ 

Tlmax = ~ °- 02 5 , linear; 

\ 

10 

2 * ’ 

7max = ~ 0.25; 


2 rt , 

^max ~ * 311 ^ 


6 . 25 * 

? ^max = 3»12 

■ 


( 11 - 8 ) 


It can be concluded that large changes occur in F(x) in the neighborhood 
of a f(x) hump. 

For large values of p, g(x) is large and sin g(x) and cos g(x) will 
oscillate with many cycles over the non-zero range of g(x). For the 
cases in Eq. (11-8) and Fig. 11-2, the approximate number of cycles in 
sin g(x) and cos g(x) for 0.4 < x < 0.8 is 1/2 for p = l /2 * and 
p = 10/2it, 2 for p = 2jt, 6 for p = 6.25*. These cycles are unequal in 
width, with the smallest width for p = 6.25* being Ax = 0.045. Since the 
interval used for the numerical integration in Fig. 11-2 was Ax ~ 0.05 9 
it is evident that the results for the p = 6 . 25 * case may not be very 
good. Either a smaller interval must be used or the numerical integration 
modified. 


To investigate this numerical integration difficulty, g(x) was 

(11-9) 


A-y Ay- 

approximated over any interval Xi - < x < Xi + — by 


where 


m = 


g(x) = g(x t ) + m(x-Xi) , 
g(x t + - g( Xl - tt) 


Ax 


( 11 - 10 ) 


If the average of t - x in Eq. (ll-2) is used for the interval At, then 
sin g(t) and cos g(t) can be integrated across the interval using the 
approximation in Eq. (11-9 )• The results are 
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i 


ti + (At/2) 


ti - (At/2) 


= - At 


sin g(t)dt = - £2ff. 8(t)l 

m Jti - 

(ti + - COS g(ti - 


(At/2) 

(At/2) 


COS g 


s(t, + f) - g(ti - f) 


= (At)Di 


( 11 - 11 ) 


f ti + (At/2) A t 

1 jH-h — 

Sin g(t t + - sin g| 

f t . _ AtVI 
r 1 2 )\ 

I too — ,riT 

J ti - (At/2) 

g( 

^ + tt) - g! 

(ti - 

At > 
2 > 

\ 


= (At)Ei . 


( 11 - 12 ) 


The term Di replaces Sj in the calculation of Ti in Eq. (11-3) and Ei 
replaces Ci in the calculation of Ui . With these terms, the results for 
the p = l/2rt, 10/2it, and 2« cases are essentially the same. However, 
for the p = 6.25n: case, the peak at x = 0.4 is about one-half as large 
and more cycles appear to be present. When twice as many intervals are 
used with the Di and Ei terms, the p - 6.25it case shows about six cycles 
in F(x) of various amplitudes and occurring in the range of the g(x) 
hump. See p = 6.25« (modified integration) case in Fig. 11-2. A similar 
cyclic result for F(x) was obtained for p = 100. 

From these results it would appear that, if F(x) shows a cyclic 
shape in a region, then there is a large deviation in the mirror in that 
region. 

In Section III for the linear case it was found a linear f(x) = C x x 
produces a constant F(x) = -2Ci . Two cases for the nonlinear equation 
were calculated by Eq. (11-4) with 


(a) g(x) = x and (b) g(x) = - ^ x (11-13) 

2it 2it 

which is equivalent to T] max = ± 1/4 at the ends. For case (a), F(x) 
varied from -8.0 at edges of the hole to -l4.8 at the ends, being 
approximately constant over the middle half of the mirror. For case 
(b), F(x) varied from + 19.6 at edges of the hole to + 4.3 at the ends, 
being approximately constant over the middle half of the mirror. 
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XII - SCALE FACTOR PROBLEMS 


Above, F(x) = I - I q is assumed known. I Q is calculated for the 
perfect mirror. I is measured on the actual mirror. If K converts I to 
scale of I 0 , then 


F(x) = 


I 

K 


lo 


( 12 - 1 ) 


It appears that K depends on the intensity of the light source and on 
the instrumentation to read I. To calibrate the perfect mirror with a 
given light of intensity, Qc. 


F(x) = 0 


K c 


/ L ipat 


( 12 - 2 ) 


For a different light intensity Q 


K = K c 
Qc 


(12-3) 


Thus, if K c can be dete rmin ed for the finished mirror ( as perfect as 
possible) using the specified instrumentation and known light Qq, then 
if the light Q can be determined, Eq. (12-3) would seem to hold, whence 


F = Sc^_ Ic 
Q, K c 


(12-4) 


A procedure to check K may be based on Eq. (12-1). If 


F(x) = Fi(x) - C 


f f F(x)dx 

J Fi(x)dx = 0 , C = h 


2(1-R) 


(12-5) 


then from Section III, a rotation of the mirror through a slope change 
of -C/2(l-R) will change I and make the new F be Fi . Thus a check K is 
(ip read after rotation), 


8l 


/ 


( 12 - 6 ) 


Jr iRdx 

k = t — 

L v* 

It may be possible to determine K directly by giving a small known 
rotation to the mirror or an equivalent displacement of the knife-edge 
and measuring the change in I at several points. Since the change in 
F(x) can be calculated (linear case) as -2C! , Ci is slope, then 

jr _(^_)ni®M. (l2— Y) 

(Al)cal 


at each point. If £E is constant, then the K may be satisfactory. 

It should be noted that an error in the origin for I shows as a 
rotation of f(x) and also an error in K shows a rotation in f(x). Thus, 
if f(x) shows a large rotation, it probably means I and K are not 
compatible, or the knife-edge is not central. 

Finally, in the Eqs. (12-2 ) and (12-6) 

f 1 82 

I I 0 dx = — Jt , no hole, 

Ji 3 1 

log* xjtEldx 

11 + X X - R 1 


for hole of radius R. 


XIII - CONCLUSIONS ML RECOMMENDATIONS 


The linear integral equation for the Foucault test of a mirror, 
solid or Cassegrain type, can be solved by several different methods. 

Since the methods of solution are approximate due to numerical integration 
and/or finite function representations, they are not unique. Although 
the methods may agree on results for some examples, they disagree on 
other examples . The inversion integral method appears to cover the 
largest range of possible problems with good results for both the solid 
and Cassegrain mirrors. The collocation method does well on many examples, 
particularly for the Cassegrain mirror, but it fails for functions with 
discontinuities. The function methods, power series and Fourier series, 
do well on some functions but have various deficiencies, particularly 
for the Cassegrain mirror. The iteration method is very poor and is not 
recommended. 


= 2tc 2 (1-R) + 2 


i 


(12-8) 
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I 


f 

1 


There appears to be several possible procedures for calibrating 
and scaling the input data for the integral equation solutions. The 
simplest procedure may be to give the knife-edge a small displacement 
and measure the change in the light intensity. Since the change can be 
calculated in the linear case, the scaling can be made by comparing the 
measured and calculated results. 

For various assumed functions for the errors on the mirror, the 
light distribution on the mirror can be calculated from the nonlinear 
equation by numerical integration. Large changes occur in the light 
intensity in the neighborhood of a local surface deviation, whether 
linear or nonlinear. However, for a large nonlinear deviation, large 
changes may occur in other regions somewhat removed from the local 
deviation. 

It is recommended that further investigation of the nonlinear 
equation be conducted, particularly with regard to 

(a) input of a nonlinear. F(x) solution into the linear 
equation and comparison of solution to original assumed 
f(x) in the nonlinear equation; 

(b) scaling and calibration in the nonlinear case; and 

(c) possible iteration solutions of the nonlinear equation. 
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APPENDIX A 


Approximation in the Sense of Least Squares 
for a System of Simultaneous Equations 

Consider the system of simultaneous equations 

N 

£ a U x i = b i 5 i = 1,2,***,M (A-l) 

d^l 

or in matrix form 

AX s B (A-2) 


If M = N, the A matrix is square and, if it is nonsingular, the system 
of equations can be solved directly. 

If M > N, the A matrix is rectangular, as there are more equations 
than unknowns. This system of equations can be solved in the sense of 
least squares in the following manner (see p. 8 of Ref. 7). 

Let xj" be a column vector that does not identically satisfy Eq. (A-l), 
but gives a residual error vector ei in 


Now 


N 

£ aijxj - bi = e 4 , i = 1,2, •••,*! 
0=1 


M M / N 

£ «? = Z £ 

i=l i=I \o=l 


\/ N 

Jy Y, a ii< x k - 



and for a minimum 


8 

CKp 



M 

/ N 

^ ^ a ip( X* " 

i=l 

\k=l 

M 

/ N 

- £ 

a ip( £ a iJ x J 

i=l 

\0=1 

M 

/ N 

2 £ 

a ip( £ a ij x J 

i=l 

\0=1 



(A-3) 


(A-4) 


(A- 5) 


p = 1,2, •••,M 
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This set of equations for a minimum of the square of the error can he 
written in matrix form as 


2[A t (AX - B)] = 0 
or 

m m ( A “6) 

A a AX = a x b 


T 

Since A A is a square matrix, symmetric, positive diagonal elements, and 
nonsingular, it can be inverted to give 


X = [A T A]" 1 A T B 


(A-7) 


as the solution of the system in the sense of least squares. 
As a simple example, consider the system 


2Xi + 3X 2 = 0 
2 X x - 4 x 2 = 1 
Xj. - x 2 =0 


Equation (A-7) gives X x = 8/45 and X 2 = - 6/45 with the residuals 


e i - “ r— 5 

45 


e * = -4F ' 


_ l4 

£3 = - — 

45 


ef + e§ + e| = - 


None of the three equations are satisfied exactly, but the X x and X 2 
values are the best values in the sense of least squares. Any other 
values of X x and X 2 will give 


e i + e I + e § > 7: 

7 

It is evident that if one or more of the equations have mistakes 
in them, these mistakes will affect all the results and residuals. The 
mistakes or errors are spread out over all the results. Such an effect 
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can occur if incompatible conditions or restrictions are present in the 
system of equations. For example, suppose the third equation in the 
simple example above is 


1000X x - X 2 = 0. 


Then, there results from Eq. (A-7) X x ~ - 4/25,000 and X 2 ~ -4/25. 
Thus the equation with the large change is still approximately satisfied, 
but the solution is completely changed. 
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APPENDIX B 


Evaluation of an Integral in Fourier Series Solution* 

The integral in Eq.. (5-11 ) 

— K mn « f* log/ , s . in Q— \ sin n0 cos m0 d0 (B-l) 

if nn J 0 6 U + cos e) 

can be written as 

2 rt 

— K ran = f log / — s - ri - - — Vsin(m+n)0 - sin (m-n)0]d0 (B-2) 

2 J n \1 + cos 0/ 


Now consider the integral 

I 2 i = f* log / - s - i - 9 V in 2i0 d0 , i = 1,2,3>**,N 
J 0 '1 + COS 0/ 

Introduce the constant l/2i and integrate by parts to get 


I 2i = 1 - cos . gi0 log sin 6 


2i 


1 + cos 0 


1” ■ ±f n m . aa 


(B-3) 


1 f 1 - COS 2i0 
2i X sin 0 


d0 


(B-4) 


Consider the difference form 


2il 2 i - (2i-2)l 2 i_ 2 = f — 210 ~ — d0 

J Q sin 0 

= - 2 Jr sin(2i-l)0d0 
o 

— [1 - cos (2i-l)it] 

2i - 1 

4 

2i - 1 

*Thahks are due to Dr. Keith Stewart son. Visiting Professor at 
The Ohio State University, Summer 19&9> and Goldsmid Professor of 
Mathematics and Joint Head of Mathematics Department, University of 
London, for this integration. 
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Since I D = 0, it follows that 


lo = 0 > 

la = - 2 , 


I4 


- — T- 

2(2) L 3 


+ 





f i 

1 k=l 



(B-6) 


Put this result in Eq. (B-6) into Eq. (B-2) and use 2i = m + n or 
m - n to get 


K nB = 


8 


(m+n)/2 

£ 


8 


|(»^)/a| 

Z-/ 


* 2 (m+n) fc-i 2k - 1 it 2 (m-n) k=l 2k -1 


(B-7) 


where m + n is even. If m = n, Eq. (B-2) gives 


m 


K„ 


-tS 


hi = n 


it m k=i 2k - 1 


(B-8) 


For m + n odd, sin (mfn)0 is even in the interval 0 to it. Since the 
log term in Eq. (B-2) is odd on the interval 0 to it, the value of the 
integral is zero for m + n odd, or 


K mn = 0 , for m + n odd . (B-9) 
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APPENDIX C 


G Matrix for Inversion Integral Method 


G MATRIX 

sm L 

0*20320 0.04411 

0.01790— 0, CL17-6.6 — 

0*02322 0*02632 

~Q • l47-01_r0. 04587- 


0.02977 

.Qjl01XZ2_ 


-0.00843 — 0*00709 
-0*00279 -0.00236 - 
ROW 2 

-0*02500 ... 0*15 8.91... 
0.01984 0.01918 

0.023.76 „ O. 02 6,79. 

-0.14677 -0.04565 - 

-0. 00831_rP. 00697. - 

—0.00273 -0.00230 - 

ROW 3 _ 

-0.00750 -0.06065 

. .0 .02255 0.02120 _ 

0.02438 0.02733 

_ - 0.1 4650- - 0.0 45.41 - 

-6.00817 -0.00685 - 
-Q. 00266 -0. 0022 5 • 
ROW 4 

— Q..J004 72 -Q.» 0.1 8 75 

0.02662 0.02402 

0.025-12 _ 0. 027.97 .. 

-0.14620 -0.04514 
-0.* 00803 -0 . 00672 
-0.00259 -0.00219 
ROW 5 

-0.00333 -0.01176 

0*0 3341_.-Q..0-28 .2 6 

0.02601 0.02872 

_ -0*14587 -0-. 04485 
-0.00787 -0.00658 

-a. 00252^-0.00212 

ROW 6 

-0*0.02-5 Q-. - Q*Q 0827 

0.04698 0.03533 

0.02709.. 0,02961 
-0.14550 -0.04452 

— 0-.0077.0_ — 0.Q0642 

-0.00244 -0.00205 

fiQyj _ 

-0.00194 -0.00618 
0.12842 0.04946 

0.02845 0.03071 

-0.14508-0.04415 
-0.00751 -0.00626 
-0.00235 -0.00198 


0.03114 
r 0.0291 a. 
'0.00602 - 
-0.00196 - 

0 ,08237 

0.01 894 
0*03156 — 
-0.02890 - 
-0. 00591.„- 
-0.00191 - 

0.16187 

0.02050. 

0.03203 

-0. 02668 ; 

-6.00580 ■ 
-0.00186 ■ 


0.02417 0. 

__CL.jQ18j05 Q. 

0.03937 0. 

r-0.X12Q76.-0, 
-0.00514 -0, 
-0.00157 -0 


02123 

.Q1869 

,05606 

.01579 

,00442 

.00113 


.0.* 03 444 
0.01906 
-0.03 9 74 
■0.02058 
•-0*00505 
-0.00153 


0,02689 
0.01953 
_ 0.05639 
-0.01563 
-0.00433. 
- 0.00110 


-0,0774.1 . 
0.02259 
0*03.258_ 
-0.02844 
-0,00568 
-0.00181 

-0.02382 
.0., 025 51- 
0.03322 
-0.02818 
-0.00555 
-0.0Q176 

-0*01489. 
0.02990 
0.03397 
-0.02789 
-0.00542 
-0. 00170 

-0.6 1042 
0.03721 
0.03488 
-0.02756 
-0. 00527 
-0.00164 


0.09606 
0,0.2031 
0.04016 
r_Q. *02-0 33_ : 

-0.00495 - 
-0,00149 - 

0,159.01 
0.02192 
._Q*Q4Q6.2_. 
-0.02017 ■ 
— 0 .90484 
-0.00145 

- 0.08881 
Q *02 4 0.7. 
0.04119 
-0,01993 
-0.00472 
— 0.Q014Q 

-0*02719. . 
0.02708 
0,04183 
-0.01967 
-0.00460 
-0.00135 

-0.01692 
0.03159 
0.04259 
-0.01938 
-0.00447 
-0. QO 130 
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0.03821 
0,02056 
0.05676 
■Q. 015.45 - 
-0.00424 - 
-0*00107 - 

0.10615 

0.02184 

0.05718 

-0.01526 

-0.00414 

-0.00104 

0.15386 

0.02349 

0.05766 

-0.01504 

-0.00404 

- 0.00101 

-0. 09766. 
0.02570 
0.05822 
-0.01481 
-0.00393 
-0.00097 

-0.02972 

0.02878 

0.05886 

-0.01455 

-0.00380 

-0.00094 


0.01951 

0.Q1968 

0.15715 

-0.01250 

-0.00380 

-0.00042 

0.02317 
0. 02039 
0.15745 
-0.01235 - 
-0,00372 - 
-0.00041 

0.02926 
0.02125 
0.15778 
-0,01219 • 
-0.00364 - 
-0.00039 

0.04146 
0.02230 
0.15815 
- 0.01201 ■ 
-0.00355 
-0.00038 

0.11462 

0.02361 

0.15858 

-0.01182 

-0.00346 

-0.00037 

0.14748 

0.02530 

0.15906 

-0.01161 

-0.00336 

-0.00036 

-0.10486 

0.02755 

0.15961 

-0.01138 

-0.00325 

-0.00034 


0.01848 

0.02112 

0.00507 

-0.01017 

-0.00326 


0.02106 

0.02174 

0.00533 

-0.01003 

-0.00319 


0.02494 

0.02247 

0.00563 

-0.00989 

-0.00312 


0.03141 

0.02334 

0.00596 

-0.00973 

-0.00304 


0.04434 

0.02441 

0.00633 

-0.00955 

-0.00296 


0.12195 

0.02574 

0.00675 

-0.00936 

-0.00287 


0.14033 

0.02746 

0.00724 

-0.00915 

-0.00277 



ROW 8 

r0a00155 — 0o00478 -0.00774 -0.01178-0.01840 -0.03170 -O 0 IIO 86 
Oo 13268 0. 13424 0.05183 0o03912 0.03340 0.03069 0.02975 

_ . ._._0 fi J)3i}L9 0*03208. 0.03^99. ...Q, 0435_0 0.05963 0.16026 0.00779 

-0.14460 -0.04373 -0.02719 -0.01905 -0.01426 -0.01112 -0.00892 

-0.007 31 -0.0 0607 -0.0051 0.. ^a^0_432__-Q . 0 0 3.67 -0.00313 -0 .00267 

-0.00226-0.00190 -0.00157 -0.00125 -0.00090 -0.00033 

_R.0W._2 . 

-0.00125 -0.00378 -0.00595 -0.00870 -0.01274 -0.01951 -0.03326 

-0^Lli2.Q__£LJL242Q-_ Qal3951 . 0^05416 Q.a0.4Lll Q. 03542 0.03295 

0.03251 0.03384 0.03737 0.04462 0.06055 0.16103 0.00844 

-0.14404 -0.04325 -0 .02 6 77_- 0a.01A6.9 -- 0^01394 -0.01084. -0.00867 
-0.00708 -0.00587 -0.00492 -0.00416 -0.00353 -0.00301 -0.00256 

-0..00217. -0.00182 -0*_Q015.0 -Q..QQ119 -Q.0QO85 -0.00031 

ROW 10 

- 9 _e.OO 1 0.2._r Q jlO Q 3 0.3 -0.00468 -0.00665. -0.00934 -0.01341 -0.02032 

-0.03446 -0.12011 0.1L653 0.14443 0.05653 0.04329 0.03776 

0.03576 0.03 619 0 .039 15 a. 0.4601 - -0^0.6.161 0.16195 0.00921 

-0.14339 —0.04269 ^0.9262,9 -0.01827 -0.01357 -0.01052 -0.00839 

-O.B-O9.6fl3_ .-Oft 00-565. -0 . 00473 -.0.00399 -0.00338 -0.00287 -0.00244 

-0.00206 -0.00173 -0.00142 -0.00113 -0.00081 -0.00030 
_RQW.JL1 

-0.00083 -0.00245 -0.00372 -0.00518 -0.00708 -0.00975 -0.01386 

-0.02088 -0.03533 -0.12359 0.10826 0.14903 0.05903 0.04577 

0.04064 0.03948 0.04152 0.04781 0.06307 0.16307 0.01013 

_ -0.14263 -0 .042Q 5 = Q. 02574 _-0. 01780 -0.01316 -0.01016 -0.00808 

-0.00656 -0.00541 -0.00451 -0.00380 -0.00321 -0.00272 -0.00231 

tP o 00 1 95 — 0 o 0 0JJ> 3. .- O. 0 Q.0 1 34 “_0_o. 001.06 -0=00076 -0.00028 

ROW 12’ 

-0.00068 -0.0 0199 - 0.002 98 -0.00 408 -0.00546 -0.00732 -0.00998 

-0.01410 -0.02120 -0.03588 -0.12637 O.IOOO’O 0.15347 0.06179 

.Q.«0_4.8_7_7__J)^04442 . 0^04485. _0,Q5_020 0..Q6487 0.16448 0,01126 

-0.14171 -0.04128 -0.02510 -0.01725 -0.01269 -0.00975 -0,00772 

-0o0Q621_-0. 00514 -0.00428 — 0,QQ35_9 -0.00303 -0.00256 -0.00217 

-0.00183 -0.00153 -0.00125 -0.00099 -0.00071 -0.00026 

ROW 13 ... 

-0.00056 -0.00160 -0.00238 -0.00322 -0.00425 -0.00557 -0.00739 

-Q.01Q02 — O.J01413 rO. 02126 -0.03611 -0.12848 0.09187 0.15790 

0.06502 0.05264 0.04983 0.05354 0.06727 0.16629 0.01267 

-0.14058 -0.04036 -0.02434 -0=01661 -0.01215 -0.00929 -0.00732 
-0.00590 -0.00484 -0.00401 -0.00336 -0.00283 -0.00239 -0,00202 

-0.00170 -0.00142 0^00116 . --0= 00092 - 0*00065 -0,00024 

ROW 14 

. -0.O-00045 -0.Q0L29 -6.00189 -Q.0Q254 -0.00330 -0.00427 -0.00554 

-0.00731 -0.00989 -0.01395 -0.02106 -0.03597 -0.12986 0.08402 

0=16255. 0 = 069.10 0.05813 0.05856 0.07064 0.16870 0.01447 

-0.13918 -0.03924 -0.02342 -0.01585 -0.01152 -0.00875 -0.00687 

-0^00^5.1 _rJI= 00450.-0=00372 -Q.Q0310 -0.00260 -0.00219 -0.00185 

-0.00155 -0.00129-0.00106 -0.00083 -0.00059 -0.00022 
_R0W 15 

-0.00036 -0.00102 -0.00149 -0.00198 -0.00255 -0.00325 -0.00416 

... -0.00537 -0.00707 -0.00957 -0.01354 -0.02056 -0.03542 -0.13044 

0.07666 0.16781 0.07474 0.06693 0.07568 0.17207 0.01689 

-Ql, 13_I37 -0.037 84: -r 0 = 02231 -0,01495 -0.01077 -0.00813 -0.00634 

-0.00506 -0.00411 -0.00339 -0.00281 -0.00236 -0.00198 -0.00166 
-0.00139 -0.00116 -0.00095 -0.00074 -0.00053 -0.00019 
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ROW 16 


- 0.00028 

- 0.00079 

- 0.00115 

- 0.00151 

- 0.00193 

- 0.00244 

- 0.00392 

- 0.00505 

- 0.00664 

- 0.00903 

- 0.01285 

- 0.01968 

- 0*13004 

Q . 07019 

0.17440 

0.08366 

0.08409 

0.17713 

- 0.13496 

- 0.03604 

- 0.02091 

- 0.01384 

- 0.00987 

- 0.00739 

- 0.00454 

- 0.00367 

- 0.00301 

- 0.00249 

- 0.00208 

- 0.00174 

- 0.00122 

- 0.00101 

- 0.00082 

- 0.00065 

- 0.00046 

- 0.00017 

ROW 17 

- 0.00021 

- 0.00059 

- 0.00085 

- 0.00112 

- 0.00142 

- 0.00177 

- 0.00278 

- Q . 0 Q 353 

- 0.00456 

- 0 . 0 Q 602 

- 0.00822 

- 0.01181 

- 0.03251 

- 0. 12832 

0.06534 

0.18405 

0.10091 

0. 18557 

, .. - 0.13158 

- 0.03364 

- 0.01912 

- 0.01246 

- 0.00877 

- 0.00650 

- 0.00394 

- 0.00316 

- 0.00258 

- 0.00213 

- 0.00177 

- 0.00147 

- 0.00103 

- 0.00085 

- 0.00069 

- 0.00054 

- 0.00038 

- 0.00014 

ROW 18 

_ - 0.00015 

- 0.00041 

— 0 . 0 Q 060 

- 0.00078 

- 0.00098 

- 0.00122 

- 0.00188 

- 0.00236 

- 0.00299 

- 0.00387 

- 0.00514 

- 0.00708 

- 0.01626 

- 0.02961 

- 0.12457 

0.06382 

0.20182 

0.20244 

- 0. 12652 

- 0.03027 

- 0.01673 

- 0.01068 

- 0.00740 

- 0.00542 

- 0.00322 

- 0.00257 

- 0.00208 

- 0.00171 

- 0.00141 

- 0.00117 

- 0.00081 

- 0.00067 

- 0.00054 

- 0.00043 

- 0.00030 

- 0.00011 

... . RJ 2 W 19 

- 0.00009 

- 0.00026 

- 0.00037 

- 0.00048 

- 0.00061 

- 0.00075 

- 0.00114 

- 0.00141 

- 0.00177 

- 0.00226 

- 0.00294 

- 0.00394 

- 0.00813 

- 0.01316 

- 0.02491 

- 0 . 11 712 

0.07071 

0.30366 

- 0.11809 

- 0.02523 

- 0.01339 

- 0.00830 

- 0.00564 

- 0.00406 

- 0.00236 

- 0.00187 

- 0.00150 

- 0.00123 

- 0.00101 

- 0.00084 

- 0 . 00057 . 

- 0.00047 

- 0 . Q 0 Q 38 

- 0.00030 

- 0.00021 

- 0.00008 

ROW 20 

- 0 ^ 0000.4 

— 0 .. Q 0012 

- 0.00018 

- 0 .00023 

- 0.00028 

- 0.00035 

- 0.00052 

- 0.00064 

- 0.00080 

- 0.00100 

- 0.00128 

- 0.00169 

- 0.00325 

- 0,00494 

- 0.00830 

- 0.01673 

- 0. 10091 

0. 16196 

- 0.10122 

- 0.01682 

- 0.00837 

- 0.00498 

- 0.00329 

- 0.00232 

- 0.00131 

- 0.00103 

- 0.00082 

- 0.00066 

- 0.00054 

- 0.00045 

- 0.00030 

- 0.00025 

- 0.00020 

- 0.00016 

- O.OOOll 

- 0.00004 

ROW 21 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

.... 0.0 

. 0,0 

0.0 

0.0 

0.0 

O.Q 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 . . . 

0.0 

0.0 

_. 0 , Q . 

0.0 

0.0 

ROW 22 

0 . 000 Q 4 

0.00011 

0.00016 

0.00020 

0.00025 

0.00030 

0.00045 

0.00054 

0.00066 

0.00082 

0.00103 

0.00131 

0.00232 

0.00329 

0.00498 

0.00837 

_ 0.01682 

0.10122 

- 0.16196 

0.10091 

0.01673 

0.00830 

0.00494 

0.00325 

0.00169 

0.00128 

0.00100 

0.00080 

0.00064 

0.00052 

0.00035 

0.00028 

0.00023 

0.00018 

0.00012 

0.00004 

..ROW 23 . 

0.00008 

0.00021 

0.00030 

0.00038 

0.00047 

0.00057 

0 * 00084 - 

0.00101 

_ o..eai23_ 

_ 0 ^ 00150 . 

_ 0.00187 

0.00236 

0.00406 

0.00564 

0.00830 

0,01339 

0.02523 

0.11809 

- 0.30366 

- 0 , 07071 . 

0.11712 

0.02491 

0.01316 

0.00813 

0.00394 

0.00294 

0.00226 

0.00177 

0.00141 

0.00114 

0.00075 

0.00061 

0.00048 

0 . Qg 037 

0.00026 

0.00009 


0.00308 

0.03433 

0.02026 

0.00572 

0.00146 


0.00222 

0.01831 

0.02533 

0.00499 

0.00123 


0.00151 

0.01030 

0.03377 

0.00412 

0.00098 


0.00092 

0.00549 

0.05066 

0.00305 

0.00069 


0.00043 

0.00229 

0.15198 

0.00172 

0.00037 


0.0 

0.0 

0.0 

0.0 

0.0 


0.00037 

0.00172 

0.15198 

0.00229 

0.00043 


0.00069 

0.00305 

0.05066 

0.00549 

0.00092 



rtOW 24 


Q* POO IX 

0.00030 

0.00043 

0.00054 

0.00067 

0.00081 

0.00098 

0.00117 

0.00141 

0.00171 

0.00208 

0.00257 

0.00322 

0.00412 

0.00542 

0.00740 

0.01068 

0.01673 

0.03027 

0.12652 

- 0.03377 

- 0.20244 

- 0.20182 

- 0.06382 

0. 12457 

0.02961 

0.01626 

0.01030 

0 *00708 

0.00514 

0.00387 

0.00299 

0.00236 

0.00188 

0.00151 

0.00122 

0.00098 

0.00078 

0.00060 

0.00041 

0.00015 


ROW 25 

0.00014 

0.00038 

0.00054 

0.00069 

0.00085 

0.00103 

0.00123 

0.00147 

Q . 001 77 

0.00213 

0.00258 

0.00316 

0.00394 

0.00499 

0.00650 

0.00877 

0.01246 

0.01912 

0.03364 

0 . 13158 

- 0.02533 

- 0.18557 

- 0.10091 

- 0.18405 

- 0.06534 

0. 12832 

0.03251 

0.01831 

0.01181 

0.00822 

0.00602 

0.00456 

0.00353 

0.00278 

0.00222 

0.00177 

0.00142 

0.00112 

0.00085 

0.00059 

0.00021 


ROW 26 

0.00017 

0.00046 

0.00065 

0.00082 

0.00101 

0.00122 

0.00146 

0.00174 

0.00208 

0.00249 

0.00301 

0.00367 

0.00454 

0.00572 

0.00739 

0.00987 

0.01384 

0.02091 

0.03604 

0.13496 

- 0.02026 

- 0.17713 

- 0.08409 

- 0.08366 

- 0.17440 

- 0.07019 

0. 13004 

0.03433 

0.01968 

0.01285 

0.00903 

0.00664 

0.00505 

0.00392 

0.00308 

0.00244 

0.00193 

0.00151 

0.00115 

0.00079 

0.00028 


ROW 27 

0.00019 

0.00053 

0.00074 

0.00095 

0.00116 

0.00139 

0.00166 

0,00198 

0.00236 

0.00281 

0.00339 

0.00411 

0.00506 

0.00634 

0.00813 

0.01077 

0.01495 

0.02231 

0.03784 

0.13737 

- 0.01689 

- 0.17207 

- 0.07568 

- 0.06693 

- 0.07474 

- 0.16781 

- 0.07666 

0.13044 

0.03542 

0.02056 

0.01354 

0.00957 

0.00707 

0.00537 

0.00416 

0.00325 

0.00255 

0.00198 

0.00149 

0.00102 

0.00036 


ROW 28 

0.00022 

Q . 00059 

0.00083 

0.00106 

0.00129 

0.00155 

0.00185 

0.00219 

0.00260 

0.00310 

0.00372 

0.00450 

0.00551 

0 . 00,687 

. 0.00875 

0.01152 

0.01585 

0.02342 

0.03924 

0.13918 

- 0.01447 

- 0.16870 

- 0.07064 

- 0.05856 

- 0.05813 

- 0. 06910 

- 0. 16255 

- 0.08402 

0.12986 

0.03597 

0.02106 

0.01395 

0. 00989 

0.00731 

0.00554 

0.00427 

0.00330 

0.00254 

0.00189 

0.00129 

0.00045 


ROW 29 

0.00024 

0.00065 

0.00092 

0.00116 

0.00142 

0.00170 

0.00202 

0.00239 

0.00283 

0.00336 

0.00401 

0.00484 

0.00590 

0.00732 

0.00929 

0.01215 

0.01661 

0.02434 

0.04036 

0. 14058 

- 0.01267 

- 0.16629 

- 0.06727 

- 0.05354 

- 0.04983 

- 0.05264 

- 0.06502 

- 0.15790 

- 0.09187 

0. 12848 

0.03611 

0.02126 

0.01413 

0.01002 

0.00739 

0.00557 

0.00425 

0.00322 

0.00238 

0.00160 

0.00056 


ROW 30 

0.00026 

0.00071 

0.00099 

0.00125 

0.00153 

0.00183 

0.00217 

0.00256 

0.00303 

0.00359 

0.00428 

0.00514 

0.00625 

0.00772 

0.00975 

0.01269 

0.01725 

_ 0. 02510 

0.04128 

0.14171 

- 0.01126 

- 0.16448 

- 0.06487 

- 0.05020 

- 0.04485 

- 0.04442 

- 0.04877 

- 0.06179 

- 0.15347 

- 0.10000 

0.12637 

0.03588 

0.02120 

0.01410 

0.00998 

0.00732 

0.00546 

0.00408 

0.00298 

0.00199 

0.00068 


ROW 31 







0.00028 

0.00076 

0.00106 

~ 0 . 0 013 4 * 

0.00163 

0.00195 

0.00231 

0.00272 

0.00321 

0.00380 

0.00451 

0.00541 

0.00656 

0.00808 

0.01016 

0.01316 

0.01780 

0.02574 

0.04205 

0.14263 

- 0.01013 

- 0. 16307 

- 0.06307 

- 0.04781 

- 0.04152 

- 0.03948 

- 0.04064 

- 0.04577 

- 0.05903 

- 0.14903 

- 0.10826 

0.12359 

0.03533 

0.02088 

0.01386 

0.00975 

0.00708 

0.00518 

0.00372 

0.00245 

0.00083 




ROW 32 

0.Q0030 0.0008L- 0.00U3_ 0.00142 0.0Q173 0.00206 0.00244 

0.00287 0.00338 0.00399 0.00473 0.00565 0.00683 0.00839 

0.Q1052 0.01357. O.OL827. 0.02629 0.Q4269 0.14339 -0.00921 

-0.16195 -0.06167 -0.04601 -0.03915 -0.03619 -0.03576 -0.03776 

—0,04329 -0*J35653_-0.L4443„-O.1L65l3 Q*J.2QLL_Q,03446_ 0, 02032 

0.01341 0.00934 0.00665 0.00468 0.00303 0.00102 

ROW. -13 „ : 

0.00031 0.00085 0.00119 0.00150 0.00182 0.00217 0.00256 

0.00301 0,.00353_. 0.00416 0.00492- 0.00587 0.00708 0.00867 

0.01084 0.01394 0.01869 0.02677 0.04325 0.14404 -0.00844 

_ -0. 1 6 103^-0.0.605 5_r,0 .04462 -0.03737 - 0 .033 84 -0.03251 -0.03295 
-0.03542 -0.04111 -0.05416 -0.13953 -0.12470 0.11590 0.03326 

0.01951 0.01274 0.00870 ... 0.00595.. 0.0Q378 0.00125 

ROW 34 

0.00033 0.00090 . 0.00125 0.00157 . 0.00190 0.00226 0.00267 

0.00313 0.00367 0.00432 0.00510 0.00607 0.00731 0.00892 

_ _Q_._01 112 0. 01 425: 0 .0 19Q5 ._0 .02719 Q ,.04273 0.1446 0 -0.0 0779 

-0.16026-0.05963 -0.04350 -0.03599 -0.03208 -0.03019 -0.02975 

-0.03069 -0._033j4Q -Q.039l2._-Q. 05183 -0^13424 -0.13268 0.11086 

0.03170 0.01840 0.01178 0.00774 0.00478 0.00155 

ROW 35 

0.00034 0.00094 0.00130 0.00164 0.00198 0.00235 0.00277 

0.00325_ .0 .003 8.0 0. Q0447 0.0052 7 0.00626 . 0. 0075 1 0.00915 

0.01138 0.01455 0.01938 0.02756 0.04415 0.14508 -0.00724 

-0.15961 -0.05885-0.04259-0.03488 -0.03071 -0.02845 -0.02746 

-0.02755-0.02878 -0.03159 -0.03721 -0.04946 -0.12842 -0.14033 

Q. 10486 .0.02972 0.01692 0.01042 0.00618 0.00194 

ROW 36 

0.00036 0.00097 _0_.QQ13.5 0,.QQ17Q„ 0.Q02Q5 0.00244 Q. 00287 

0.00336 0.00393 0.00460 0.00542 0.00642 0.00770 0.00936 

0.01161 0.01481 .0.01967 0.Q2789 0.04452 0.14550 -0.00675 

-0.15906 -0.05822 -0.04183 -0.03397 -0.02961 -0.02709 -0.02574 

-0,02530 -0.02570 -0,02708 -0.02990 -0.03533 -0.04698 -0.12195 

-0.14748 0.09766 0.02719 0.01489 0.00827 0.00250 

ROW 3-7 •- - ... - - 

0.00037 0.00101 0.00140 0.00176 0.00212 0.00252 0.00296 

0.00346 0.00404 0.00472 0.00555 0.Q0658 0.00787 0.00955 

0.01182 0.01504 0.01993 0.02818 0.04485 0.14587 -0.00633 

-0.15858 -0.05766 -0.04119 -0.03322 -0.02872 -0.02601 -0.02441 
-0.02361 -0.02349 -0.02407 -0.02551 -0.02826 -0.03341 -0.04434 

r_0.*_l 146Z. -0,1533.6 0,-08881... 0,02382- 0.01176 0.Q0333 

ROW 38 

__ 0.00038 -0.00104 0.00145 Q.0.018L ... 0.Q.0219 0.00259 0.00304 

0.00355 0.00414 0.00484 0.00568 0.00672 0.00803 0.00973 

0.01201 0.01526 0.02017 Q. 02844 0.04514 0.14620-0.00596 

-0.15815-0.05718 -0.04063 -0.03258 -0.02797 -0.02512 -0.02334 
_ -0, 0223 0 -0 .021 84 -0 .021 S2__.-..Q ,Q225.9__- Q . 024Q2 -0.02662 -0.03141 

-0.04146 -0.10615 -0.15901 0.07741 0.01875 0.00472 

ROW 39 ■ , - 

0.00039 0.00107 0.00149 0.00186 0.00225 0.00266 0.00312 

. 0,00364 0.00424 .0.00495 0.Q058Q. 0,00685 0.00817 0.00989 

6.01219 0.01545 0.02038 0.02868 0.04541 0.14650 -0.00563 

-0.15778 -0.05676 -0.04016 -0.03203 -0.02733 -0.02438 -0.02247 

-0.02125 -0.02056 -0.02031 -0.02050 -0.02120 -0.02255 -0.02494 

-0.02926 -0.03821 -0.096Q6 -0.16187 0.06065 0.00750 

95 


ROW 40 

— - 0* 0 0041 QjjQQl.lQ_ 0*00153 0.00191 0.0.0230. 0. 00273 . 

0.00372 0.00433 0.00505 0.00591 0.00697 0.00831 

0.01235 0.01563 0.02058 0. 02890 0. 04 5 6 5. 0*14677. 

—0.15745 -0.05639 -0.03974 -0.03156 -0.02679 -0.02376 

-0.02039 -0. 0.1953 -0.0 1906 -0. 0L894 -0*01918 -Q. 01984 

-0.02317 -0.02689 -0.03444 -0.08337 -0.15897 0.02500 

ROW. 41. _ ___ : 

0.00042 0.00113 0.00157 0.00196 0.00236 0.00279 

0.00380 0.00442 0.00514 0.00602 0 .00709 Q.Q0843_ 

0.01250 0.01579 0.02076 0.02910 0.04587 0.14701 

-0. 15715 -0. .056 06 _^0 . 03937_ -0 * 03 114 -0*02632. -0.02322 
-0.01968 -0.01869 -0.01805 -0.01772 -0.01766 -0.01790 - 
. — 0. 01951 -0.02123 -0.02417 -0.02977 -0.04411 -0 . 2 032 0 


Q. 00319. 
0.01003 
-J0..0Q533L 
-0.02174 
-0«021Q6_ 


0.00326 

0.01017 

-0.00507 

- 0.02112 

-0.01848 
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APPENDIX D 


Computer Program for Inversion Integral 
Method (Cassegrain Mirror) 

The computer program listed here is a Fortran program written for 
the IBM SYSTEM 360-75 . The program solves the integral equation for the 
Foucault Test using the inversion integral method of solution for the 
case of a mirror with a central hole. All operation are performed in 
double precision (approximately 1 6 decimal places),, 

The main program reads in M and MODE under a (215) format . M is 
the number of evenly spaced points on the mirror and must be even. The 
end points and points on the edge of the hole are among these M points, 
MODE determines which method of numerical integration is to be used 
(MODE = 1, trapezoid rule; MODE = 2, Simpson’s rule; MODE = 3? 5-point 
rule). Also in the main program R, the radius of the central hole, is 
read in under a (E 12,8) format. Next the subroutine GMA.T is called. 
GMAT generates the G matrix which is stored in G(l, J). The subroutine 
GMAT calls subroutine AMA.T, which calculates the A(I, J) matrix 
(Eq. 10-36). MAT in turn calls the subroutine HMAT, which calculates 
the numerical integration matrix H(l). H(l) is used in A(l, J) and 
A(I, J) is used in calculating G(l, J), Eqs. (10-35), (10-33), (10-37), 
and (10-40). Control then returns to the main program where NRUN, the 
number of cases to be solved, is read in under an (15) format. Next the 
subroutine CAFF is called. CAPF calculates the light intensity differ-* 
ence F(x) from a functional representation at the M points on the Mirror 
(Eqs. 10-42, 10-44, 10-50, 10-51) o The M values of F(x) are stored in 
FBIG(l) in the main program and F(l) in subroutine CAPF. Subroutine 
DGMPRD then matrix multiplies F into FBIG to return FLIT, Eq. (10-39) « 
The matrix FLIT contains M values of the surface error f (x) calculated 
at the same points as F(x). DGMFRD is an external IBM subroutine con- 
tained in the Scientific-Subroutine Package. If the ’hump case" for 
f(x) is to be run, as in case 4 of this listing, A and B, the location 
and size of the hump, must be given values in CAPF. If F(x) is to be 
inputed data rather then functionally represented, then CAPF may be 
replaced by a READ statement. 
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c 


IMPLICIT RJ=AL *8 (A- -1 ,0-7.) 

0 I M»=fJS I , ' 1 “‘ A ( *2 » 42 ) , H ( 4? ) , G ( 42 , 42 ) , F 8 10(42) » FL I T ( 4? ) 
jM^IN'/FRS ip.-i I '.'X F G R A L 'FTHOO FOR wn LF C'\ SF*** 

Ul R I T F ( 6 , 1 0 9 ) 

109 FORMAT ( / , 3. OX , 1 ***Ii'VF RS ION INTEGRAL OFTHOD FOR HOi_ F CASE***') 
RPAO( *5,101 ) ?J,«or)F 

IfU FORMAT (2 15) 

REAO( 5 , 102 ) R 
10? FORMAT (3E12. 8) 

ORITE(6,107) »' , MMO F , R 

107 FORMAT (///,10X, 'MO. OF POIMTS = ' , 1 5 , / , 10X , * > O'lE = » » I *5 , 1 3 X , » R 
*,E16.8) 

0 A L L G H A T ( k , K nn F , G , R ) 

MRITE( 6, 104) 

104 FORMAT ( /// , 10X , 1 G MATRIX ' ) 

no z 1 = 1 

2 1 ‘RITE (6,103) I , (G( I ,j ) , J=1 ,M ) 

3.03 FORMAT ( 1H , 5X , ' RO'-i i , I 3/ ( 3 .OX , 5E 16 .« ) ) 

REA 0(5, 101) : 'R|JM 
MRITE(6, 13.0) MR l IN 

110 FORMAT ( 10X, 'MO. nF CASES RUN = * * 1 5 ) 

00 3 1=3. »NRUN 

CALL CAPF (M. ,FRIG ,R , I ) 

CALL OG'iP'->0(G, FRIG, FLIT, 1) 

!| RITE(6,105) 

103 FORMAT (///, 10X, 'LITTLF F') 

3 R I T F ( 6,106) ( F L I T ( J ) , J =3. , M) 

106 FORMAT ( 3. H , ( 10X , 5F1 6 .8 ) ) 

STOP 

F- %, 0 
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SUBROUTINE GMAT ( M , M 00 E ,G,R ) 

IMPLICIT REAL*8 ( A-H»G-Z ) 

DIME MS ION G (42 ,42 ) »H (42 ) ,0 (42 ) »GT (42 ) , A( 42, 42 ) , X (42 ) 

P 1=3.1415926535897900 
P I2=P I *P I 

CALL AMAT(M,MnDF,A,H,R,X) 

DO 1 J=1,M 

1 D( J)= PI2+(A(J,J))**2 

MH=H/2 

NH1 = '1H+1 . 

U1=M-1 

NHM1=MH-1 

MH2=MH+2 

D ( 1 ) =28 .000 

D ( M) =28 «0D0 

D ( 1 1H ) = 1 1 .0 DO 

D( MH1 ) =11 .000 

DO 5 1 = 1, M 

DO 2 J = 1 , M 

I F ( J . FO . I ) GO TO 2 

G ( I , J ) =X ( J > *A ( I , J ) / ( X ( I ) *0 ( J ) ) 

2 CONT INIJE 

G ( 1,1) =— A ( I ,1 )/Q ( I ) 

5 CONTINUE 

G (1,1) =0.20 6000 
G(M,li) =-0.20 6000 
G ( MH ,MH ) =-0 .31 80 DO 
G(MH1,MHI) =0.318000 
G( 2,1) =-0.028300 
G(M1 , : -i) =0.02 8300 
G ( i i Hi i 1 , M H ) = 0 . 0 43 2 DO 
G(i-.H2 ,‘-Hl )=-n. 0*3200 
G ( r-iMi'il ,’iHl ) =-0 o 00 ? 1 DO 
G(l-iH2 ,: ; H ) = 0. 003100 
on 9 1=1, R 
H ( I)=G(1H,I ) 

9 GT( T ) =G ( TH1 , I ) 

DO 8 1 = 1, (i 
DO 8 J=l,’1 

8 G ( I , J ) =G ( I , J ) - ( ( 1 . ODO-R / ( X ( I ) ) )*H( J) + ( 1 . OO0+ 1/ ( 

* ) ) / 2 .one 

RETURN 

END 


XU))) *GT ( J 
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S t IB ROUT I ME A MA T ( M , MODE , A , H , R , X ) 

IMPLICIT RFAL*8 (A-H,n-Z) 

n I MENS I ON A ( M , M ) , H ( 1 ) , A T ( 42 » 42 ) , 8 ( 4 2 , 42 ) , X ( 42 ) 
MH=w/2 

CALL HMA T ( H H » M 00 F , H ) 

!'■' 1 = M- 1 
M 2 
ImH 1 = MH+1 
MM2 = .)H+2 
i 'H3=Mh+3 

Mh , j'!2 = MH-2 

P = 2*( 1.000-R)/DPL0AT (N>2 ) 
no i j=i,mh 

1 x(j)=-i.ono+p*'(j-i) 

no 2 J= 'HI ,i- 

2 X( J)=R+P*(J-MH1 ) 
on 3 1 = 1 ,»•' 

no 3 J = 1,M 

IF ( J »EO. I ) GO TO 3 
A ( I , J ) = ( P*H ( J ) ) / ( X ( J ) -X ( I ) ) 

3 CONTINUE 

00 4 1 = 2, "‘l 

IP1=I+1 

I - '1=I-1 

A { I , I'll )=-H{ INI )-0.500*H{ I ) 

A { I » I P 1 ) =H ( I P 1 )+0.500*H( I ) 

A no= o.o no 

DO 5 J=1,M 

I F ( I .FO. J ) GO TO 3 
A 00= A 00-r A ( I , J ) 

5 CONTI i '• IE 

4 o(I,I) =-A00 

A ( i ;H , .-’Hi ) = ( p*H C : H1 ) ) / ( X ( ‘-’H I ) -X ( MH ) ) 

A ( i >H 1 ,i‘K>=(P*H(*'H> )/ ( X ( H ) -X (MHJ. ) ) 

A ( i H, *>H «1 )= (P#H(MHM1 ) ) / (X ( ■ • H * * 1 ) -X C-!H ) ) 

A ( ,- IJ l , 5 H2 ) = ( P*H ( -«H2 ) ) / ( X ( MH2 ) -X ( KH l ) ) 

A ( 1 , 1 ) = 1 . 0 no 
A ( =1 .000 

"RITE (6,: 1.0? ) 

102 F 0 R ' A T ( / / / , 1 0 X , * A V, ATR I X ’ ) 

00 10 1 = 1,-' 

10 '-‘RITE (6,101) I,(A(I,J),J = 1, :; ) 

101 F n R'iAT ( 1H ,5X, , 13 , / ( 10X ,5E 16 .8 ) ) 

RETURN 

END 
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S l IS R OUT I E K "A T ( 00 c , |-; ) 

I ! i P L I C IT R E A L * 8 ( A. : - H , 0 - 7. ) 

1 1 If i E.^'S 10 i ! I { * l ) 
i-' l = i-l 
■?. = \ -2 

r;« i to (i ,2 ,3 ) » ’nn!f 

1 nn io T = 2 ,‘"-i 
10 H{ I) =1.000 

|-'( 1 ) = 0 a 5 1)0 
H(fi)=O.Hi)0 

fin to 7 

2 on 11 1 = 2 , 1 ''!,? 

U K( I ) =*-.()00/3 .000 
OH 12 1=3,02,2 
12 l*!( I) =2. 000/3. 000 
H ( 1 ) = 1 .000/3.000 
M ( } = l.ono/3. non 
GO TO 7 

3 on 16 1 = 2,; '1,2 

16 l-i( I) =32.0 00/^-5 .0 00 
00 17 1=3,02,4 

17 l-l( I) =12. 000/45. 000 
;i4=;;-4 

on 1.8 1 = 5 ,04,4 

18 H ( I ) = 1 4 . 0 00 / A 5 . 0 00 
H ( 1 ) = 7.000/45.000 
H( 01=7.000/45.000 

7 Cni'ITIHJE 
on 4 l=i,i 
i-;i=i-:-I 

4 H ( i ll ) =H ( I ) 

,v .2 = M+ - -t 

'•RITE ( 6 , 101 ) 

101 FORMAT ( / // , 10X , *H MATRIX 1 ) 

WRITE ( 6, 102 ) (1-1(1 ) ,1=1,02) 

102 FORMAT (1H , ( 10X , 5E16. 8 ) ) 
RETURN 

EMD 
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SI!:* ROUTINE C A °P ( =•' , c > K , J ) 

I i PI. IC IT RF/>L*3 (A-M,n-Z) 
i> Iii*=NS IOi* F (I-') 
f.'? = v /2 

p=2.ono*(i .ono-R)/ (?'- 2 .onn) 

X=-l .nno-p 
A = 0 o 6 5 11 0 
'•=0.17500 
"RITE ( 6 , 3.03 ) -\ , « 

103 FPR'-'AT (// »10X. ' a = ' , F16.H,/, 10X, ' ■* = ' ? E 15 » I ' ) 
on 52 1 = 1 j- ' 

X=X+P 

on TO ( 1,2,3 ,4,5,0 ,7 ) , J 

1 F ( I ) = -30 .OOfs 
GO TO 52 

2 F< I ) =-35 ,onn*x 
GO TO 52 

3 F { I ) =- 3 5 . 0 DO *X *X 
GO TO 52 

2= ( X-/' ) / o 

IP (7** 2 .GI-. l.or.r.nj) r;n To 53 

IF (2**2. .IT. 1.00001 ./‘-r’-n. 7. ** 2 ,GT. . 99990 ) on In 55 
r ( I ) = — 2 o 0 — 7_**3 + n.O .0/3 . 0 ) *7- ( ( 1 . 0-7** 2 ) ** 2 ) *‘>U IG ( IV-.P.S (((1.0- 
10+2) )*( (1 ,0+X )/ (1 ,0-X) )*( (X-K) / ( X +R ) ) ) ) 

GO TO 52 

5 3 Ft I)=2.0/ ( 3 .0*2 ) +2 .0/ ( 3 . n*z**3 ) -( ( 2** 2- 1 . 0 ) **? ) * ( 2. 0/2 + 2.0/ ( 
1*3) + 0I.OG (OAFS ( ( 7-1 .0 )/ ( 2 + 1.0 ) ) ) ) 

Gil TO 52 

55 F( I)=HSIG '( (4. 0/3.0 )0) ,2) 

GO TO 52 
5 F ( I ) = 1 4 . 0 O0 / X 
GO TO 52 

o F ( I ) = OS I GO ( 3 5 . 0 no , x ) 

GO TO 52 

7 F( I ) =l*.ono/ (X*X) 

52 IF (I .FO. f+2) X = R-P 

!'RITF(6, 1.01 ) (F ( I ) , 1 = 1 , '• ) 

101 FOR -'AT ( // » 10X » 'RIG F ' , / , ( 60 1 6 . 3 ) ) 

RFTORO 
/ FHi) 


7 ) / ( 1 
3 . ■ ) * 7 


102 
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